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1.0  SUMMARY  OF  RESEARCH  OBJECTIVES 


During  the  period  from  January  1  through  December  31,  1983  we 
continued  our  research  In  the  following  two  areas: 

1.  Continued  experimental  studies  of  soft  x-ray  population 
inversions  in  recombining,  laser  produced  plasmas  with  an 
emphasis  on  the  development  of  a  practical  soft  x-ray 
amplifier. 

2.  Analysis  of  the  performance  which  could  realistically  be 
expected  to  be  obtained  from  multilayer  structures  suitable 
for  use  as  normal  incidence,  soft  x-ray  cavity  mirrors. 

Our  activies  during  this  period  have  led  to  several  significant 
achievements: 

1.  Operation  of  our  frequency  tripled  laser  system  in  line- 
focus,  step  target  geometry. 

2.  Observation  of  significant  population  Inversion  in  line 
focus  experiments  on  the  Balmer  <*  transition  in  hydrogenic 
flourine. 

3.  Successful  numerical  modelling  of  both  point-  and  line-focus 
experiments. 

4.  Completion  of  the  analysis  of  possible  soft  x-ray  multilayer 
mirror  performance. 

Two  theses  were  successfully  completed  during  this  period  and  a 
summary  presentation  of  the  line-focus  experimental  results  was  given 
at  the  conference  on  Laser  Techniques  for  Extreme  Ultraviolet  Spectros¬ 
copy.  Since  a  complete  description  of  the  work  is  contained  in  the 
theses,  they  are  attached  and  constitute  the  main  body  of  this  report. 
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ABSTRACT 


The  feasability  of  an  XUV  amplifier  based  on  the  recombination  of  a 
laser-produced  plasma  was  investigated.  Experiments  were  carried  out 
with  a  Nd+^-glass  laser  generating  sub-nanosecond  pulses  of  either  IR 
light  (1.054  ym,  up  to  100  J)  or  frequency- tripled  UV  light  (0.351  pm, 
up  to  60  J) .  The  plasma  formed  at  the  surface  of  a  flat  target  was 
cooled  at  an  early  stage  of  its  expansion  by  inserting  a  metallic  foil 
near  the  target.  The  cooling  of  a  dense  plasma  is  known  to  result  in  a 
recombination  cascade  dominated  by  collisions,  which  tend  to  populate 
high- lying  quantum  states  of  the  recombining  ion. 

With  point-focus  illumination  (diam.  <100  pm),  population 

inversions  were  observed  in  the  helium-like  ions  Al+11  and  Mg+1°, 

showing  possibility  of  gain  on  the  transition  3  D-4  F  at  129.7  A  and 

154.4  A  respectively.  We  established  that  the  inversions  were  caused 

neither  by  resonance  mechanisms  such  as  photopumping  or  charge  transfer, 

nor  by  channelling  of  the  plasma  flow.  We  defined,  after  a  parametric 

study,  the  optimum  target  design  (slit-foil  of  an  arbitrary  material, 

200  pm  wide  and  200  pm  from  the  original  plasma  formation) ,  and 

14  2 

illumination  conditions  (UV  laser  of  intensity  =5x10  W/cm  ) .  The 
short  wavelength  laser  was  shown  to  couple  its  energy  to  the  target  more 
efficiently. 

With  line-focus  illumination  (100  pmx2  mm,  using  specially 
manufactured  cylindrical  optics),  inversion  was  observed  in  the 
hydrogenic  ion  F  between  levels  n«2  and  n«3.  The  corresponding 
transition,  Balmer  a  at  80.9  A,  was  spatially  resolved  over'vl  mm  in  the 


direction  of  expansion,  with  single-shot  recording. 

Nuaerical  simulations,  using  existing  laser-fusion  code  and  rate 
equations,  were  in  good  agreement  with  the  data  when  we  included 
enhanced  recombination  (via  cooling)  and  plasma  opacity.  The  gain 

0 

region  travelled  with  the  velocity  of  the  expanding  plasma  (0.3-1  10 

2  -1 

cm/sec).  The  gain  could  be  high  locally  in  time  and  space,  10  cm 
over <  100  psec  in  Al+11  (129.7  k)  and  5  cm-1  over  ^200  psec  in  F4®  (80. 
k) .  However,  the  plasma  length  (200  pm  and  2.5  mm  respectively,  in  the 
region  of  highest  gain)  did  not  permit  direct  gain  measurement  from 
time-integrated  recordings  of  the  amplified  spontaneous  emission.  We 
discuss  how  that  goal  could  be  reached  in  the  future  with  upgraded  UV 
pump  laser  (factor  5)  and  time-resolved  diagnostic  capability  ($100 
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Chapter  1 
Introduction 


.•H 


The  interest  in  building  an  X-ray  laser  comes  from  the  potential 
applications  of  a  coherent  source  of  X-rays  in  fields  like  chemistry, 
biology  and  crystallography  (LLL  1974-75).  Sources  of  high  brightness, 
emitting  subnanosecond  pulses  of  energetic  photons  would  make  possible 

studies  combining  the  spatial  resolution  of  X-rays  with  a  temporal 
resolution  similar  to  what  is  now  available  in  the  visible  and  near  UV. 
Some  applications,  like  diffraction  experiments  in  biology,  require  only 
a  point  source.  But  spatial  coherence  would  provide  directionality  to 
the  emission,  an  important  factor  in  a  spectral  region  where  imaging  is 
limited  to  small  solid  angle,  grazing  incidence  optics. 

Both  components  of  such  a  laser,  the  amplifying  medium  and  the 
resonator,  present  us  with  technological  problems  yet  to  be  solved. 

Among  the  ideas  for  producing  the  amplifying  medium,  a  promising  one  is 
the  rapid  generation  of  a  highly  ionized  plasma  by  focusing  on  a  solid 
target  the  intense  pulse  from  a  subsidiary  laser.  Such  an  approach, 
which  benefits  directly  from  the  advances  in  laser  development  for 
^nertial  confinement  fusion  (ICF)  research,  has  been  pursued  at  the 
Laboratory  for  Laser  Energetics  (LLE)  for  several  years  (Forsyth  et  al 
1976,  Bhagavatula  fc  Yaakobi  1978).  The  present  thesis  describes  the 
progress  made  in  that  direction,  including  experimental  results  and 
numerical  simulations.  A  separate  study  of  the  resonator  is  also  under 
way  at  LLE  (Rosenbluth  1982).  Available  pump  laser  energies  limit 
current  research  to  the  XUV  region  ('vlOO  k) ,  not  much  farther  than  the 


region  already  accesible  by  frequency  multiplication  (Reintjes,  et  al, 
1981)  or  mixing  (Egger,  et  al,  1982)  of  existing  lasers.  A  direct 
approach  is  attractive,  however,  because  it  is  often  scalable  to  shorter 
wavelengths. 

1-1.  Inversion  methods. 

There  are  two  basic  difficulties  in  producing  an  amplifying  medium 
operating  in  the  X-ray  region  (LLL  1974-75).  The  first  is  the 
photoelectr ic  effect  which  results  in  the  absorption  of  X-rays  except  in 
hot,  highly  ionized  plasmas.  The  second  is  the  very  short  lifetime,  T, 
of  allowed  X-ray  transitions,  which  places  severe  requirements  on  the 
pumping  mechanism  sustaining  the  population  inversion: 

T(sec)  =  (8Tr2r#c)-1(gu/g1)A2/flu  =  1.5  x  10“ 1  6  (g^/g^)  [A  (A)  ]  2/f^ 

where  rQ  is  the  radius  of  the  electron  (2.82  10~5  k) ,  c  is  the  velocity 
18 

of  light  (3  10  k/ sec) ,  X  is  the  wavelength  of  the  transition,  g^  and 

g^  are  the  statistical  weights  of  the  upper  and  lower  levels 
respectively,  and  f  is  the  absorption  oscillator  strength  (Griem  1964, 
p.45).  f^  is  nearly  constant  along  isoelectronic  sequences,  and 
approaches  unity  for  intense  lines  (Wiese,  et  al,  1966-69). 

Some  authors  (Lacour  et  al  1976,  Duguay  1976)  have  proposed  pumping 
schemes  fast  enough  to  be  effective  in  a  time  shorter  than  T  for  the 
production  of  ultra-short  bursts  of  coherent  X-rays.  In  such  devices, 
an  unambiguous  diagnostic  would  be  difficult  due  to  the  low  photon 
yield.  Recent  progress  in  conventional  laser  technology  has  made 
possible  a  steady-state  type  of  pumping,  while  providing  conditions  such 
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N(X),  PARTICLE  DENSITY  DISTRIBUTION 


N  (x) 


LASER  PULSE 


1_.  Recombination  to  Hydrogenic  Plasma  (i.e.  one  bound  electron). 


A.  Collisional  Processes: 

three-body  recombination 
and  its  inverse 

collisional  excitation 
and  de-excitation 

ion-ion  collisions 

B.  Radiative  Processes; 

radiative  recombination 
and  photoionization 

spontaneous  emission 

absorption 

(optically  thick  plasma) 


,  z+  - 
(N  +e  )  +  e 


-►nIZ.^  +  +  e  (accelerated) 
—  (p) 


_(z-l)+  -  superelastic^  „(z-l)+ 

N  .  ,  +  e  .  - - ■»  N  +  e 

(p)  inelastic  (q)  (q<p) 


NZ+  +  e" 


■*N;z:1)+  +  hv 

—  (p) 


N 


(z-l)+ 

(P) 

(z-l)  + 


■>N(q) 


(Z-l)+ 


n;  +  hu 

(q) 


-*■  N 


+  hv 

(z-1) + 

(p) 


(p  >  q) 


—•  Recombination  to  lower  stages  of  ionization. 

It  involves  processes  analogous  to  those  listed  above,  and  also 


dielectronic  recombination 
and  its  inverse  (Auger) 

charge  exchange  (here  one 
example  among  many) 


(z-l)+  - 

N  ,  .  +  e 

P 


-*■  N 


2  N 


(z-l)+ 

(p) 


-►  N 


(z-2) + 

(P1'P2) 
(z-2) + 
(P1»P2) 


+  N 


(pl'p2  >  p) 
z+ 


3_.  Free  electrons. 

bremsstrahlung 
and  its  inverse 


e 


*  e  (slowed)  +  hv 


Table  1-1  :  Plasma  Recombination  Mechanisms, 


as  to  prevent  severe  absorption  effects.  When  a  high-power  laser  is 

focused  on  a  target,  it  creates  a  plasma  at  the  surface  (Figure  1-1). 

The  plasma  remains  in  a  state  of  non- equilibrium  for  a  time  of  the  order 

of  the  laser  pulse  duration,  xL»  which  can  be  much  longer  than  x. 

During  that  time,  x^,  recombinations  in  the  plasma  take  place  through 

mechanisms,  both  collisional  and  radiative,  which  populate  the  various 

levels  of  the  lower  stages  of  ionization  at  different  rates  (Table  1-1) . 

In  collisional  processes  an  ion  can  interact  with  another  ion  (with 

possibilities  of  charge  exchange),  one  electron  (inelastic  collision, 

dielectronic  recombination)  or  two  electrons  (three-body  recombination) . 

Radiative  processes  can  involve  one  ion  (spontaneous  emission, 

stimulated  emission  and  absorption) ,  one  ion  and  one  electron  (radiative 

recombination)  or  one  electron  (bremsstrahlung) .  Some  of  the  ideas 

suggested  to  make  this  plasma  behave  as  an  X-ray  amplifier  have  led  to 

experimental  works:  charge  exchange  (Scully,  et  al,  1973),  resonant 

photoexcitation  (Norton  &  Peacock  1975) ,  dielectronic  recombination 

(Jaegle,  et  al,  1981),  and  recombination  with  cooling  (Gudzenko  & 

Shelepin  1965) ,  the  latter  being  the  basis  of  our  own  research. 

(i)  Charge  exchange.  The  idea  is  to  utilize  electron  pick-up  in  an 

excited  state  of  the  ionized  species.  Dixon  et  al  (1977-78)  have  taken 

advantage  of  the  large  cross-section  for  the  transfer  of  an  electron 
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from  atcxnic  Carbon  to  C  '  resulting  in  the  selective  production  of 
excited  ions  [C**'+*)  .  They  focus  a  Nd+^-glass  laser  (5  J,  16  nsec)  on 
a  carbon  slab.  A  background  gas  (Helium,  1  torr)  is  used  to  allow  the 
formation  of  neutral  atoms  and  to  slow  down  the  reacting  ions,  a  low 
relative  velocity  being  required.  Population  inversion  occurs  between 


the  levels  n-5,4,3  in  both  and  C-^  about  2  cm  from  the  slab  but  the 
low  ion  density  ('ilO  cm  )  results  in  such  a  weak  emission,  and  weak 
gain,  that  many  shots  are  necessary  for  measurement  and  no  diagnostic  is 
possible  on  the  potential  amplification  (350-760  A  region) .  In  order  to 
achieve  higher  densities,  Seely  &  McKnight  (1977)  have  proposed  to  use 
Argon  instead  of  Helium,  the  charge  transfer  taking  place  directly 

I  £  ^2 

between  the  ions  C  and  Ar  within  1  mm  of  the  target.  Amplification 
is  expected  at  182  A  on  C+^(3-2).  No  experimental  follow-up  on  this 
scheme  has,  so  far,  been  reported. 

(ii)  Resonant  photoexcitation.  In  a  dual  plasma  the  radiation  emitted 

by  one  species  ("source  ion")  pumps  the  other  species  ("lasing  ion")  to 

an  excited  level  n^2  (n,  principal  quantum  number),  resulting  in  an 

inversion  between  the  levels  n  and  n-1.  The  process  is  effective  if 

there  is  resonant  absorption  of  the  pump  radiation,  as  can  be  done  by  an 

appropriate  choice  of  species.  Hagelstein  (1981)  has  proposed  to  pump 

Ne+^  "by  preferential  photoionization  of  2s  electrons  from  [atomic 

neon] ,  resulting  in  a  2s- 2p  inversion  at  27  eV  [460  A]  in  Ne+1.  The 

target  design  implementing  the  scheme  is  based  on  thin  film  flashlamp 

and  filter  systems  which  emit  strong  X-ray  radiation  near  500  eV  [25  A] 

14 

when  irradiated  with  short  pulse  (90  psec)  1.06  ym  radiation  at  10 

2  -1 
W/cm  .  (...)  Gains  in  excess  of  50  cm  are  calculated"  (ibid.). 

Experiments,  to  be  carried  out  at  LLE  with  the  OMEGA  laser  system 

(Bunkenberg,  et  al,  1981),  are  in  the  planning  stage. 

(iii)  Dielectronic  recombination.  Ions  whose  configuration  includes 
more  than  one  excited  electron  often  correspond  to  energy  levels  above 
the  ionization  limit,  that  are  called  "autoionizing"  because  they  can 


decay  spontaneously  by  ejecting  an  electron  instead  of  by  emitting  a 
photon.  The  populations  of  these  levels  are  in  equilibrium  with  the 
thermal  continuum  of  free  electrons  even  when  the  ionic  distribution 
itself  departs  from  equilibrium.  Thus  the  reverse  process  of  electron 
capture  (dielectronic  recombination)  followed  by  collisional  transfer  to 
levels  below  the  ionization  limit  may  enhance  strongly  the  populations 
of  some  levels  of  specific  ions,  and  be  a  source  of  population 
inversion.  This  is  the  interpretation  given  by  Jaegle  et  al  (1981)  for 
their  own  early  observation  of  gain  at  117.41  A  on  the  transition 
2p54d(3P1)-2p6{1So)  in  Al  IV  (1971,  1974).  A  gain  coefficient  of  10 
cm  1  is  reported  in  a  plasma  of  electron  density  lO^-lO21  cm-3  formed 
by  a  Nd+3-glass  laser  (20  J,  40  nsec). 

(iv)  Recombination  with  cooling.  Collisional  and  radiative  processes 
tend  to  favor  high-  and  low-lying  quantum  levels  respectively.  In  a 
dense  plasma,  collisions  become  important  and  an  inversion  can  be 
achieved  during  the  cascade  through  the  excited  levels  of  the 
recombining  ions.  Rapid  cooling  is  required  in  order  to  "freeze"  the 
plasma  in  a  much  higher  ionization  state  than  predicted  by  the 
equilibrium  relations.  "Cold"  electrons  recombine  mainly  into  the 
higher  levels,  while  the  lower  levels  are  still  relatively  empty 
(Gudzenko  t  Shelepin  1965) .  Population  inversions  have  been  observed  in 
particular  in  Carbon  VI  during  the  adiabatic  expansion  of  a  cylindrical 
fiber  (Jacoby  et  al  1981) ,  and  in  Aluminum  XII  through  external  heat 
loss  (Bhagavatula  6  Yaakobi  1978).  In  the  former,  a  direct  measurement 
of  gain  x  length  product  up  to  5  is  reported  on  the  3-2  (Balmera  ) 
transition  at  182  A.  In  the  latter,  the  cooling  is  so  effective  that 
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line-intensity  inversions  are  observed  and  a  gain  of  10  cm-1  is 

3  3 

predicted  on  the  manifold  4  F-3  D  at  129.7  %. 

1-2.  Thesis  outline. 

The  plasma  recombination  concept  is  the  most  promising  of  those 
listed  above:  (i)  it  has  already  been  demonstrated  in  the  visible 
(Silfvast  et  al  1979-b,  1980) ,  (ii)  several  independent  experiments  at 
various  laboratories  have  led  to  population  inversion  in  the  XUV  region, 
(iii)  there  is  potential  for  significant  gain  because  inversion  can 
occur  at  high  density.  Our  thesis  being  based  on  that  concept,  in 
Chapter  2  we  outline  the  theoretical  treatment  of  recombining  plasmas 
developed  by  Bates  (1962-a  and  -b) ,  McWhirter  (1963)  et  al.  We  also 
review  several  diagnostic  methods  relevant  to  this  thesis. 

Chapter  3  gives  some  details  about  the  instrumentation  necessary 
for  an  experimental  study:  high-power  laser,  focusing  optics,  diagnostic 
instruments.  For  the  most  part,  we  used  existing  hardware,  with  some 
modifications  to  fit  our  requirements. 

Our  study,  described  in  Chapter  4,  tried  to  address  three  issues: 
the  choice  of  an  XUV  transition,  the  target  design  to  achieve  a 
population  inversion,  and  the  design  of  an  amplifier  for  direct  gain 
measurement.  It  is  desirable  to  choose  a  transition  from  a  highly 
stripped  ion,  to  limit  the  number  of  recombination  channels.  Also,  the 
higher  the  atomic  number,  the  shorter  the  wavelength  of  the  transition 
is,  but  an  upper  limit  is  set  by  the  available  pump  laser.  Experiments 
were  carried  out  with  a  Nd+^-glass  laser  generating  sub-nanosecond 


pulses  of  either  IR  light  (1.054  ym,  up  to  100  J)  or  frequency-tripled 


OV  light  (0.351  ym,  up  to  60  J) . 

With  point-focus  illumination  (diam.  £100  ym) ,  population 
inversions  were  observed  in  the  helium-like  ions  Al+11  and  Mg+1°, 
showing  possibility  of  gain  on  the  transition  33D-43F  at  129.7  X  and 
154.4  A  respectively.  We  established  that  the  inversions  were  caused 
neither  by  resonance  mechanisms  such  as  photopumping  or  charge  transfer, 
nor  by  channelling  of  the  plasma  flow.  We  defined,  after  a  parametric 
study,  the  optimum  target  design  (slit-foil  of  an  arbitrary  material, 

200  ym  wide  and  200  ym  from  the  original  plasma  formation) ,  and 

m  a  *y 

illumination  conditions  (UV  laser  of  intensity  =5x10  W/cm  } .  The 
short  wavelength  laser  was  shown  to  couple  its  energy  to  the  target  more 
efficiently. 

Experiments  were  also  done  with  line-focus  illumination  (100  ymx2 
mm,  using  specially  manufactured  cylindrical  optics) ,  to  study  the 
design  of  an  actual  XUV  amplifier.  Inversion  was  observed  in  the 

48 

hydrogenic  ion  F  between  levels  n*2  and  n«3.  The  corresponding 
transition,  Balmer  a  at  80.9  A,  was  spatially  resolved  over  ^1  mm  in  the 
direction  of  expansion,  with  single-shot  recording. 

We  made  numerical  simulations  of  our  results  with  a  1-dimensional 
laser  fusion  code,  LILAC  (Goldman  1973,  Delettrez  1978-82).  Chapter  5 
shows  how  the  hydrodynamic  motion  of  the  plasma  (LILAC) ,  the  ionic  rate 
equations  (McWhirter  &  Hearn  1963) ,  and  our  own  radiation  transport  and 
target  modeling  were  combined.  Experimental  data  could  be  reproduced 
when  we  included  enhanced  recombination  (via  cooling)  and  plasma 
opacity.  The  gain  region  travelled  with  the  velocity  of  the  expanding 

fi 

plasaa  (0. 3-1.0  10  cm/sec) .  The  gain  could  be  high  locally  in  time  and 


space,  10*  cm  over  <100  psec  in  A1  (129.7  %)  and  5  cm"  over  ^200 
psec  in  P  (80.9  A).  However,  the  plasma  length  (200  ym  and  2.5  mm 
respectively,  in  the  region  of  highest  gain)  did  not  permit  direct  gain 
measurement  from  time-integrated  recordings  of  the  amplified  spontaneous 
emission.  That  goal  could  be  reached  in  the  future  with  upgraded  UV 
pump  laser  (factor  5)  and  time-resolved  diagnostic  capability  (  <100 
psec) . 

More  details  about  the  instrumentation  and  data  reduction 
procedures  can  be  found  in  three  appendices.  They  cover  our  photometry 
&  film  calibration  (Appendix  A),  the  experimental  wavelength 
determination  of  lines  recorded  with  crystal  and  grating  spectrographs 
(B) ,  and  the  design  &  fabrication  of  the  cylindrical  optics  to  form  a 
line-shaped  plasma  for  an  XUV  amplifier  (C) .  The  work  described  in  the 
appendices,  while  an  important  part  of  our  experimental  study,  can  also 
be  useful  to  spectroscopists  in  other  applications.  Several  users' 
codes  have  been  written  to  make  that  work  easily  available. 


Chapter  2 


Recombination  of  a  Highly-Ionized  Plasma 


There  are  a  number  of  proposals  to  form  an  X-ray  amplifying  medium 
in  a  recombining  plasma  {see  Chap.l).  For  a  more  quantitative 
evaluation  of  their  feasability,  a  comprehensive  calculation  of  the 
population  densities  in  the  plasma  is  required.  The  standard  model  is 
that  of  McWhirter  &  Hearn  (1963).  Modifications  and  refinements  have 
been  suggested  for  some  coefficients  (Whitney  et  al  1974-1980,  Weisheit 
1975,  Landshoff  &  Perez  1976,  Salzmann  &  Wendin  1978),  but  the  overall 
procedure  is  widely  accepted.  In  this  chapter  we  outline  their  model, 
which  is  the  basis  of  the  atomic  rate  equations  routine  included  in  our 
numerical  simulation  (Chapter  5) .  We  also  review  some  spectroscopic 
techniques  to  measure  the  temperature  and  density  in  a  highly-ionized 
plasma. 


2-1.  Rate  Equations. 

McWhirter  &  Hearn  have  calculated  the  instantaneous  population 

densities,  P(n),  of  the  excited  levels  (n,  principal  quantum  number)  of 

hydroqenic  ions  in  a  plasma.  The  steps  in  their  procedure  are  (a)  to 

get  an  expression  for  the  rate  coefficient  of  each  collisional  or 

radiative  process,  (b)  to  write  the  rate  equations  for  the  P(n),  (c)  to 

make  reasonable  assumptions  in  order  to  simplify  the  equations  and  solve 

them,  for  the  electron  temperature  (T  )  and  density  (N  )  prevailing  in 

e  e 

the  plasma.  The  collisional  and  radiative  processes  listed  in  Table  1-1 
are  all  included  except  those  involving  the  absorption  of  photons: 


opacity  is  not  considered,  at  least  not  at  this  stage.  One  looks  for  a 

local,  "quasi-static"  solution  for  a  given  set  of  conditions  (T  ,n  ) . 

e  e 

The  assumptions  that  are  made  are  the  following: 

(i)  The  free  electrons  have  a  Maxwellian  distribution. 

(ii)  Steady-state  approximation.  "The  time  in  which  the  population 
densities  of  the  excited  levels  come  into  equilibrium  with  a  particular 
population  density  of  the  ground  level,  free  electrons  and  bare  nuclei 
is  so  short  that  it  is  sufficient  to  express  these  equilibrium 
population  densities  as  functions  of  the  ground  level  population 
density"  (ibid.).  Setting  the  time  derivatives  dp(n)/dt=0  (n>l)  leads 
to  a  set  of  simultaneous  equations, 

P(n)  =  C  +  C.  P(l)  (n  >  1)  (2-1) 

o  1 

where  Cq  and  are  time-independent  parameters  which  are  functions  of 

n,  N  ,  T  and  the  atomic  number  Z. 
e  e 

(iii)  The  number  of  simultaneous  equations  (2-1)  can  be  kept  finite 
because  the  field  from  the  surrounding  plasma  produces  a  lowering  of  the 
ionization  limit:  beyond  a  certain  value  of  n,  the  bound  levels  lose 
their  discrete  nature  and  merge  into  the  continuum  of  free  electrons. 
Also,  for  levels  of  high  ,i,  collisions  dominate  and  the  levels  are  in 
Saha-Boltzmann  equilibrium  (which  is  simply  the  result  of  Boltzmann 
statistics  applied  to  the  system  [ion  +  electron],  see  for  example  Griem 
1964) .  That  approximation  is  made  for  n>20  "because  this  is  the  limit 
of  the  tables  of  oscillator  strength.  (...)  For  many  conditions  this  is 
much  higher  than  necessary"  (McWhirter  6  Hearn  1963).  Eq.2-1  can  be 
written  in  terms  of  the  Saha-Boltzmann  population  densities  P  (n) 


defined  for  the  local  value  of  (T  ,N  )s 

e  e 

P (n) /P£ (n)  *  rQ(n)  +  ^ (n) [P (1) /PE (1) ]  (2-2) 

One  still  has  to  solve  the  rate  equation  for  the  ground-state.  It 
can  be  written  as  follows: 

dP(l)/dt  =  NgtaP(N+Z)  -SP(1)1  (2-3) 

where  P  (N42)  is  the  population  density  of  bare  nuclei.  In  that 

"collisional-radiative"  model,  the  recombination  coefficient  a  includes 

three-body  and  radiative  recombination,  collisional  deactivation  and 

cascading,  while  the  ionization  coefficient  5  includes  collisional 

excitation,  and  ionization  and  photoionization,  a  and  S  are  functions 

of  T  ,  N  and  Z. 
e  e 

A  further  generalization  is  possible  if  one  defines  "reduced" 
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temperature  and  density  ne*N@/Z  .  The  quantities  a/ Z  and  z  S 

are  functions  of  and  ne  only.  Hence  the  important  parameters  to 

consider  when  studying  a  highly  ionized  plasma  are  (6  ,n  ) ,  not  (T  ,N  ) 

e  e  e  e 

McWhirter  &  Hearn  have  generated  on  a  computer,  in  the  range  0.35 
8  —3  18  —3 

eV<0<22  eV  and  10  cm  <n<10  cm  ,  the  coefficients  r  (n)  and  r,  (n) 
e  'e  o  1 

(see  Eq.2-2),  and  the  ratio  P  (1)/P  (1)  where  P  (1)  corresponds  to  a 

S  E  S 

true  steady-state  (ground  level  in  equilibrium). 

In  the  extreme  cases  of  high  and  low  density  ne»  the  recombi nation 
is  controlled  by  collisional  and  radiative  processes  respectively. 
Collisions  help  fill  the  upper  excited  states  of  the  hydrogenic  ion, 
while  radiative  decay  rapidly  populates  the  ground-state  and  lower 
excited  states.  The  filling  of  the  upper  levels  through  collisions  is 


also  influenced  by  temperature.  Let  us  consider  the  three-body 
recombination  and  its  reverse  process  (Table  1-1) .  In  a  collision 
between  a  free  electron  and  an  hydrogenic  ion,  ionization  is  more  likely 
to  occur  if  the  bound  electron  is  already  in  an  upper  level.  This  is 
especially  true  when  the  incident  electron  has  a  low  kinetic  energy, 
i.e.  low  temperature  6@.  Seen  from  the  opposite  perspective,  the 
probability  of  a  collision  between  one  ion  and  two  free  electrons  is 
higher  at  lower  velocity.  Most  interactions  will  be  at  a  distance, 
producing  weak  energy  perturbations,  and  recombi nation  will  occur  in  one 
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Plasma  Soectrosco 


Plasma  temperature  and  density  may  be  inferred  from  the 
spectroscopic  recording  of  transitions  in  1-  and  2-electron  ions  (Griem 
1964,  1974).  Among  the  lines  that  are  useful  for  such  diagnostic, 
and/or  are  present  in  many  of  our  experimental  results  (Chap. 4),  we  will 
mention  the  following: 

(i)  Resonance  lines.  They  are  allowed  transitions  to  the  ground 
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state,  ls-np  and  Is  (  S)-lsnp(  P)  (n£2) ,  in  the  hydrogenic  and 

helium-like  ions,  respectively.  The  series  ls-np  is  also  called  Lyman 

series  (Lyman  a  for  n=2,  0,Y/<S»e*etc. . .  for  n=3,4,5,6,etc. . . ) .  When  the 

term  "resonance  line"  is  used  alone,  it  refers  to  the  first  line  of  the 

2 

series,  i.e.  ls-2p  or  Is  -ls2p,  often  the  latter  when  there  is  no 
ambiguity. 

(ii)  Intercombination  lines.  They  are  forbidden  transitions  to  the 

2  1  3 

ground  state,  such  as  the  triplet-singlet  lines  Is  (  S)-lsnp(  P)  and 
2  1  3 

Is  (  S)-lsns(  S)  in  the  helium-like  ion.  In  our  case  the  label  will 
always  refer  to  n=2. 

(iii)  Satellite  lines.  They  are  transitions  corresponding  to  a 

doubly-excited  level,  as  can  be  produced  by  dielectronic  recombination. 

For  example,  lines  of  the  Lyman  series  (ls-np)  have  helium-like 

satellites  lsms-npms  and  lsmp-npmp  (m£2) .  Similarly,  lines  of  the 

2 

helium-like  resonance  series  (Is  -lsnp)  have  lithium-like  satellites 
2  2 

Is  ms-lsnpms  and  Is  mp-lsnpmp.  The  most  intense  satellites,  and  those 


resulting  in  a  spectral  shift  large  enough  to  be  resolved,  correspond  to 

2  2 

m«2,  especially  ls2p-2p  (labeled  s)  and  Is  2p-lsnp2p  (labeled  s^  unless 
specified  otherwise).  There  are  also  higher  order  satellites,  where  the 


same  basic  transition  takes  place  in  an  ion  whose  electronic  shells  are 

2 

partially  filled.  The  last  satellite  of  the  resonance  line  Is  -ls2p  is 
the  so-called  line,  when  the  atom  is  neutral  except  for  one  vacancy 
created  in  its  K-shell. 

2-2-a .  Temperature  Measurement. 

One  may  distinguish  between  methods  which  require  some  knowledge  of 

the  electron  density  N  and  those  which  do  not  (Woodall,  et  al,  1976). 

e 

One  finds  in  the  first  category  (i)  the  ratio  of  two  resonance  lines 
2 

(ls-np)/(ls  -lsnp) ,  especially  with  n*2,  (ii)  the  ratio  of  two  lines 

belonging  to  the  same  resonance  series.  The  second  category  includes 

(i)  the  ratio  of  a  resonance  line  and  its  nearby  satellites,  which  is 

proportional  to  Tfi,  (ii)  the  slope  of  the  free-bound  continuum  beyond 

the  series  limit,  which  is  proportional  to  exp(-hv/kTe) .  A  serious 

problem  in  dense  plasmas  is  that  of  self-absorption.  Methods  involving 

lines  with  high  oscillator  strength,  hence  high  opacity  (such  as  ls-2p 
2 

and  Is  -ls2p) ,  become  very  inaccurate.  Examples  will  be  given  later 
(Figs. 5-4  6  5-5). 

2-2-b.  Density  Measurement. 

Available  methods  include  (i)  the  ratio  of  the  helium-like 
resonance  to  intercombination  lines,  (ii)  Stark  broadening.  The  first 
method  relies  on  the  fact  that,  in  a  certain  range  of  densities,  the 
lifetimes  of  the  resonance  (S)  and  metastable  (3p)  states  are 
respectively  shorter  and  longer  than  the  time  constant  for  collisions 


between  the  two  states.  That  time  constant  ,  hence  the  line  ratio,  is 


density  dependent  (Griem  1975) .  In  the  second  method,  the  spectral 
broadening  induced  by  fields  from  surrounding  ions  makes  the  discrete 
spectrum  merge  with  the  free-bound  continuum.  The  number  of  the  last 
discernible  spectral  line  in  the  resonance  series  depends  on  the 
broadening,  which  is  density  dependent  (Inglis  &  Teller  1939). 


Instrumentation 


Experiments  were  carried  out  with  several  laser  configurations,  in 
terms  of  wavelength,  beam  diameter,  output  energy  and  power.  Optics  to 
form  a  point  focus  or  a  line  focus  were  used,  sometimes  with  the 
addition  of  a  blast-shield.  Diagnostic  instruments  included  a  variety 
of  crystal  spectrographs,  grating  spectrographs  and  pin-hole  cameras. 
This  chapter  describes  that  instrumentation  and  the  context  in  which  it 
was  used. 

3-1.  Laser . 

i 

Figure  3-1-a  shows  the  Glass  Development  Laser  ,  GDL  (Seka  et  al 
1980),  a  one-beam  prototype  of  the  24-beam  OMEGA  laser  built  for 
Inertial  Confinement  Fusion  (ICF)  research:  "[its]  basic  design 
characteristics  include  Nd-doped  phosphate  glass,  rod  amplifiers  up  to 
90  mm  diameter,  the  propagation  of  a  circularly  polarized  beam  [to 
reduce  the  effective  non-linear  index  of  refraction  of  the  glass] , 
extensive  spatial  filtering  and  imaging,  and  the  use  of  large  aperture 
Pockels  cells"  (ibid.).  Depending  on  the  type  of  oscillator  used,  the 
output  at  A«1.054  ym  can  be  a  "short",  high-power  pulse  (50  psec  FWHM, 
up  to  500  GW)  or  a  "long",  energetic  pulse  (500  psec,  up  to  100  Joules; 

1  nsec,  up  to  120  Joules) .  These  numbers  correspond  not  to  peak 
performance  but  to  what  the  system  can  deliver  on  a  regular  basis.  Also 
available  to  us  for  some  experiments  was  an  up-graded  version  of  GDL 
where  a  series  of  prototype  high-power  amplifier  modules  called  "active 
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mirrors"  (Brown  et  al  1981)  were  added  at  the  end  of  the  beamline.  They 
boosted  the  output  of  the  laser  above  200  Joules  (700  psec  FWHM) ,  with  a 
beam  diameter  of  150  mm. 

In  the  past  two  years  GDL  has  been  converted  to  short  wavelength 
(0.351  ym)  by  frequency  doubling  and  tripling  the  1.054  ym  beam  in  a 
pair  of  type  II  KDP  crystals  (Seka  et  al  1981) .  The  90  mm  GDL  output 
beam  is  expanded  to  130  mm  to  nearly  fill  the  crystals.  The  overall 
conversion  efficiency  to  third  harmonic  is  close  to  60%  in  long  pulse 
operation  (500  psec) .  The  layout  is  shown  in  Figure  3-1-b.  A  number  of 
diagnostics,  mainly  calorimeters  and  streak  cameras,  provide  a  complete 
characterization  of  the  laser  pulse. 

3-2.  Optics. 

The  GDL  output  beam,  at  1.054  ym  or  at  0.351  ym,  is  linearly 
polarized.  A  series  of  coated  mirrors,  optimized  for  maximum 
reflection,  transport  the  beam  toward  the  focusing  optics  and  the 
experimental  vacuum  chamber.  When  the  direct  IR  output  is  used,  a 
quarter-wave  plate  (mica)  is  added  at  the  end  of  the  beamline  so  that 
light  reflected  by  the  target  cannot  propagate  back  and  damage  the  rod 
amplifiers.  That  risk  does  not  exist  in  the  UV  configuration  because 
the  frequency  conversion  is  not  reversible. 

The  optics  used  to  focus  the  beam  on  target  and  the  flat  window 
serving  as  interface  with  the  vacuum  are  made  of  BK-7  glass  (1.054  ym 
light)  or  fused  silica  (0.351  ym  light),  with  anti-reflection  coatings 
on  each  face.  In  all  point-focus  studies  a  single-element  aspheric  lens 
was  used.  In  early  IR  experiments  the  lens  had  a  focal  length  of  30  cm, 
a  working  f-number  3.5  and  was  located  inside  the  vacuum  tank.  In  order 


to  go  through  its  clear  aperture,  the  90  mm  output  beam  from  GDL  was 
slightly  demagnified,  to  85  ram,  with  an  inverse  telescope.  In  later  IR 
experiments  the  lens  had  a  focal  length  of  60  cm  and,  depending  on 
whether  or  not  the  active  mirror  booster  amplifier  was  used,  a  working 
f-number  4  or  6.  That  lens,  located  outside  the  tank,  and  the 
vacuum-interface  window  were  similar  to  the  focusing  optics  used  in  the 
OMEGA  laser.  Designed  as  a  unit,  they  are  corrected  for  spherical 
aberration  and  coma.  In  UV  experiments  the  lens,  located  outside  the 
tank,  had  a  focal  length  of  167  cm  and  a  working  f-number  12.7. 

Line-focus  studies  were  also  carried  out,  using  cylindrical  optics. 


In  IR  experiments  a  single-element  double-aspheric  lens  produced  a  1  mm 

line  in  the  first  astigmatic  image  plane  (focal  length  f*30  cm, 

f  -f  *0.35  cm,  working  f-number  3.5).  In  UV  experiments  two  weakly 
y  x 

diverging  cylindrical  correctors  were  added  before  the  main  focusing 
lens  used  in  point-focus  studies.  Appendix  C  gives  details  about  the 
design  of  the  cylindrical  optics  and  describes  how  the  insertion  of  a 
blast  window  in  the  beam  converging  on  target  affects  the  distribution 
of  intensity  in  the  focal  plane.  The  blast  window  is  intended  to 
protect  from  target  blow-off  the  focusing  optics  or  the  vacuum  interface 
window,  whichever  is  directly  exposed. 

3-3.  Diagnostics. 

We  relied  mainly  on  soft  X-ray  spectroscopy  to  study  the  plasma 
created  by  the  laser  beam  at  the  surface  of  the  target.  Two  types  of 
transitions  were  of  interest  to  us.  The  transitions  which  end  in  the 
ground  state  of  the  Hydrogen-like  or  Helium-like  ion  give  a  means  to 


measure  the  populations  of  the  excited  states  of  these  ions  and, 
combined  with  some  physical  assumptions,  the  electron  temperature  and 
density  of  the  plasma.  For  the  materials  we  studied,  these  transitions 
range  from  ^  1.5  A  (Nickel,  atomic  number  Z=28)  to  above  20  A  (Oxygen, 
Z*8) .  They  were  recorded  with  crystal  spectrographs.  We  were  also 
interested  in  the  transitions  ending  in  excited  states  of  these  same 
ions.  They  correspond  to  the  lines  which  are  candidates  for 
amplification  by  stimulated  emission.  In  suitable  materials  (from 
Aluminum,  Z*13,  to  Fluorine,  Z*9)  their  wavelengths  are  in  the  100  A 
range.  They  were  recorded  with  a  grating  spectrograph.  In  some 
experiments  a  pin-hole  camera  was  added  to  monitor  the  size  and  aspect 

ratio  of  the  plasma.  This  was  especially  relevant  when  a  line-focus  was 
used.  The  experimental  set-up,  shown  in  Figure  3-2,  took  place  in  a 
large  (17"high  x25"diam)  vacuum  chamber  at  /v  10-5  torr. 


3-3-a.  Crystal  spectrograph. 

X-rays  emitted  by  the  plasma  are  diffracted  according  to  Bragg's 
law  A*2D.sin@ , where  0  is  the  angle  of  incidence  on  the  diffracting 
planes  of  the  crystal  and  D  is  the  interplanar  spacing.  The  spectrum 
is  recorded  on  film.  Appendices  A  and  B  give  details  about  the 
photometry  of  the  spectrograph  and  its  dispersion  characteristics, 
respectively.  As  shown  in  Figure  3-2,  it  is  usually  equipped  with  a 
set  of  filters  and  a  slit  parallel  to  the  direction  of  dispersion  of  the 
crystal.  The  filters  are  thin  foils  which,  in  addition  to  providing  a 
spectral  selection,  serve  as  an  entrance  window,  preventing  light  leaks 
inside  the  instrument.  Figure  3-3  gives  the  transmission  curves  of  the 


Figure  3-2  :  Experimental  Set-up  for  XUV  Inversion  Studies. 
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filters  we  used  most  frequently:  Beryllium,  Aluminum,  and  Saran 
^C2H2C^2^x*  sPat*a*  resolution  is  provided  (i)  along  the  normal  to  the 
target  by  the  slit  (typical  width  ^50  ym) ,  (ii)  in  the  transverse 
direction  by  the  crystal  itself  due  to  its  small  acceptance  angle  at  any 
specific  wavelength.  The  slit  can  be  either  a  pair  of  semi-infinite 
razor-blades  or  two  stripes  (or  rods).  The  latter  case  is  outlined  in 
Figure  3-4,  where  the  laser  is  focussed  on  a  flat  target  through  an 
opening  formed  by  two  foils  (see  Chap. 4).  The  hydrogenic  and 
helium-like  resonance  series  of  the  target  material  are  recorded  on  the 
film.  Five  regions  can  be  identified,  from  left  to  right:  (I) 
integrated  spectrum  of  the  whole  plasma,  (II)  shadowgram,  (III) 
spatially  resolved  spectrum  produced  by  the  slit,  (IV)  inverse 
shadowgram,  (V)  integrated  spectrum  of  part  of  the  plasma.  Radiation 
from  the  dense  region  of  the  plasma  near  the  surface  of  the  target  does 
not  reach  the  film  on  the  right  side  of  the  slit  because  it  is  blocked 
by  the  edge  of  the  target.  Weaker  emission  from  the  under-dense  region 
is  not  swamped  and  its  features  can  be  observed  in  (IV) .  Also  the 
spectrum  on  the  right  side  (V)  corresponds  to  a  spatial  integration 
limited  to  the  expanding  plasma.  Table  3-1  summarizes  the  main  features 
of  the  spectrograph  for  each  of  the  materials  studied.  Crystals  with 
long  2D-spacings,  TAP  being  one  example,  have  a  rather  large  acceptance 
angle  roughly  proportional  to  the  wavelength.  The  spatial  resolution 
provided  by  the  crystal  thus  decreases  when  spectra  from  low-z 
materials  are  recorded.  For  a  typical  distance  of  8  cm  between  the 
source  and  the  crystal,  the  spatial  resolving  element  exceeds  100  ym  in 


an  oxygen  plasma. 


Crystal-Spectrograph  Character i st i cs . 


Grating  spectrograph.  (Figure  3-5) 


The  grating,  laid  on  a  concave  surface  of  radius  r,  not  only 
disperses  the  incident  XUV  radiation  but  also  forms  a  one-dimensional 
image  of  the  entrance  slit,  on  a  curved  film  (radius  r/2) .  All  three 
elements,  the  slit,  the  grating,  and  the  film,  are  on,  or  tangent  to, 
the  so-called  Rowland  circle  of  radius  r/2  (Sawyer  1963) .  That 
arrangement  insures  that  the  imaging  condition  is  satisfied,  independent 
of  wavelength,  because  the  angle  of  incidence  a  is  constant  across  the 
grating.  An  additional  slit,  in  a  direction  perpendicular  to  the 
entrance  slit  and  located  anywhere  between  the  source  and  the  film,  can 
provide  one-dimensional  spatial  resolution  of  the  plasma  emission,  with 
arbitrary  magnification.  The  design  of  the  spectrograph  must  meet  the 
following  requirements: 

a/  Spectral  range.  It  must  cover  a  band  of  at  least  100  A  centered 
around  the  line  of  interest,  for  a  complete  identification  of  the 
various  ionic  species  in  the  plasma  (see  Appen.B) . 

b/  Resolution.  The  thermal  Doppler  broadening  gives  an  upper  limit 
for  the  realizable  resolution  of  two  adjacent  lines: 


(  A/AA)  *  410  [M,  ,/T.  .  J 

D  (amu)  i(kev)J 


(3-1) 


where  N  is  the  ion  mass  and  T^  the  ion  temperature.  Since  typically  T^ 
ttlOO  ev  in  our  experiments,  it  proved  very  difficult  to  study  some 
promising  laser-transition  candidates  such  as  Al+3,3l<33D-43F)  at  129.7  X 
(AAd=.03  X)  which  overlaps  with  the  strong  Al+3(2p-3d)  line  at  129.73  A. 
Other  transitions  are  more  isolated.  The  nearest  neighbors  to  the 
Mg+l0(33D-43F)  line  at  154.4  A  are  Mg44 (2s  2p5-2s  2p4  3s)  at  152.5  A  and 
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44  a 

possibly  0  (2p-5d)  at  153.95  A,  if  the  target  used  is  free  of  Iron  and 

Chromium  impurities.  The  Balmer  lines  in  F  are  very  easily  resolved. 
Since  we  are  only  trying  to  record  particular  transitions,  without 
studying  their  line-profile,  a  resolving  power  51^1000  is  appropriate 
(AA  M).l  K) . 

c/  Etendue.  It  is  difficult  to  set  absolute  specifications  because  in 
the  EUV  region  there  is  no  reliable  calibration  for  the  film  we  used, 
Kodak  101-01,  and  no  data  for  the  diffraction  efficiency  of  commercial 
gratings.  However  one  can  compute  the  relative  photometric  efficiency 
of  a  particular  geometry  (see  Appen.A) . 

d/  Dimensions.  The  optimal  geometry  of  the  spectrograph  for  a  given 
wavelength  is  determined  by  the  choice  of  grating  (rao: us  of  curvature, 
spatial  period,  blaze  angle).  Additional  constraints  are  introduced  by 
the  vacuum  environment  (inner  diameter  of  the  tank  and  of  its  ports) . 

A  code  (filename:  DESIGN)  has  been  written  to  select  among  the 
commercially  available  gratings  those  which  are  compatible  with  both  the 
optical  and  mechanical  requirements.  For  any  grating  the  code  computes 
(i)  the  angles  of  incidence  and  diffraction  which  satisfy  the  blaze 
condition  at  the  wavelength  of  interest,  (ii)  the  spectral  dispersion 
and  resolution,  (iii)  the  relative  photometric  efficiency,  (iv)  the 
dimensions  of  the  spectrograph  and  how  they  fit  in  the  tank.  For  our 
purpose  the  most  suitable  gratings  were  the  Bausch  &  Lomb  models 
35-5200-400,  710  and  810,  of  radius  of  curvature  400.7  mm.  Their  number 
of  grooves/mm  and  blaze  angle  are  (600,  1°43'),  (1200,  3°33')  and  (2400, 
6°54'),  respectively.  They  offer  various  compromises  between  proximity 
to  the  blaze  condition  (model  400)  and  good  resolution  (model  810) . 


Models  710  and  810  were  ordered  with  a  gold  coating  to  improve  their 
efficiency.  They  also  had  the  advantage  (after  reduction  of  their  blank 
size)  of  being  ready  for  use  in  a  rocket  spectrograph  lent  to  us  by  the 
Naval  Research  Laboratory,  which  contains  a  film  cassette  of  the 
appropriate  curvature. 

The  entrance  slit  was  set  to  provide  an  angle  of  incidence  a  -  5°. 
Its  width  was  about  15  vim  and  its  distance  to  the  source  4-5  cm.  no 
window  was  required:  an  aperture  stop  and  proper  shielding  were  enough 
to  prevent  light  leaks.  Eastman-Kodak ,  the  manufacturer  of  101-01  film, 
suggests  that  a  distilled  water  presoak  could  minimize  infectious 
development  effects.  We  found  that  it  did  produce  a  more  uniform,  less 
granular  film  but  the  spectral  lines  were  more  fuzzy.  We  decided 
against  presoaking,  to  avoid  a  loss  of  spectral  resolution. 

3-3-c.  Pin-hole  camera. 

A  pin-hole  camera  was  used  in  a  limited  number  of  experiments  to 
measure  the  size  and  aspect  ratio  of  the  X-ray  emitting  plasma.  It  was 
made  of  a  <t>  4  pm  hole  located  d=l  cm  from  the  source  and  f=10  cm  from 
the  film  plane  (Kodak  RAR-2491) .  A  25  ym  Beryllium  filter  limited  the 
observation  to  the  region  \  <  10  &  (see  Fig. 3-3).  To  avoid  vignetting, 
the  observation  was  made  at  45°  from  the  normal  to  the  target  (which  was 
also  the  laser  axis),  providing  a  2-dimensional  image  of  magnification 
10  x  The  condition  for  optimal  resolution  (Woodall  et  al  1976) , 

$  ( (d+f )/Xdf ]  -2,  is  satisfied  for  \  =5  The  corresponding  resolution 
element  is  A  *5/6  (1+  d/f)<J><  .  This  was  typically  the  spectral 

region  we  were  interested  in  (see  Chapter  4). 


Experimental  Results 

Recombination  lasers  in  a  laser-vaporized  target  material  have  been 
made  to  work  in  the  visible  part  of  the  spectrum  (Silfvast  et  al  1979-b, 
1980).  In  order  to  extend  their  range  to  EUV  and  X-ray  wavelengths,  one 
needs  either  to  consider  higher  stages  of  ionization  or  to  move  to 
higher  atomic  numbers  Z  along  an  isoelectronic  sequence.  How  far  one 
can  go  depends  on  the  available  pump,  in  our  case  a  laser  at  1.054  pm  or 
0.351  pm.  This  chapter  presents  results  on  the  interaction  between 
laser  and  matter,  which  were  the  basis  for  our  choice  of  target 
materials  and  XUV  transitions  to  investigate.  A  study  of  the  target 
design  was  carried  out  (i)  to  optimize  a  population  inversion  in  the 
XUV  domain,  (ii)  to  understand  the  mechanism (s)  leading  to  that 
inversion.  Experiments  were  done  on  the  development  of  an  actual 
device,  which  could  lead  to  a  meaningful  gain  measurement  from 
recordings  of  the  amplified  spontaneous  emission  (ASE) . 

4-1.  Laser-matter  interaction. 

A  laser  of  intensity  I  is  focused  on  a  target  of  atomic  number  Z. 

13  15  2 

Depending  on  the  value  of  I  (10  to  3x10  W/cm  ) ,  a  dense,  highly 
ionized  plasma  can  be  created  at  the  surface  of  targets  of  z-number  as 
high  as  ^30.  All  electrons  are  stripped  except  those  of  the  inner 
shell  (K-shell) ,  which  have  higher  ionization  barriers.  The  optical 
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electric  field  (>10  V/cm)  not  only  exceeds  the  breakdown  value  for 


solid  materials  but  produces  strong  acceleration  of  the  plasma 
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electrons,  resulting  in  an  electron  temperature  in  the  kilovolt 

range.  The  ionization  stage  one  reaches  has  an  ionization  potential  X 

of  the  order  of  3  to  5  times  that  temperature  (Whitney  &  Davis  1974) , 

provided  the  plasma  is  heated  long  enough  for  ionization  equilibrium  to 

n 

be  established.  X  is  proportional  to  z  while  T  has  a  relatively  slow 

e 

dependence  on  I  (T^  ^  I^3) .  Hence  there  is  a  very  steep  scaling  for 
the  intensity  necessary  to  reach  the  K-shell  of  targets  of  increasing 
Z-number . 

The  laser  light  actually  propagates  through  the  plasma  only  up  to 

the  point  where  the  natural  plasma  oscillation  frequency, 

w  *(4ire2N  /m)Vi,  is  equal  to  the  laser  frequency  w  (N  ,  electron 
p  €  L  6 

density,  e  and  m,  charge  and  mass  of  the  electron  respectively,  all  in 

CGS  units)  (Hora  1981).  This  happens  when  N  (cm-3)  =1021/ (XT  (ym) ) 2, 

called  "critical  density"  (Ne)c>  The  strong  dependence  of  (Ng)c  on  the 

wavelength  X  of  the  incident  light  shows  that  one  may  expect 

significant  variations  in  the  coupling  of  different  lasers  with  a 

plasma,  hence  the  interest  in  doing  experiments  at  two  different 

wavelengths  (1.054  ym  and  0.351  ym) . 

At  short  wavelength  the  laser  radiation  penetrates  to  higher 
22  3 

density  (10  electrons/cm  ) ,  where  the  plasma  is  more  collisional  and 
where  collisional  absorption  ,  i.e.  inverse  bremsstrahlung ,  is  more 
effective.  At  1  ym  there  is  also  resonant  absorption  (resonant  coupling 

1  4  2 

between  light  and  plasma  waves)  when  the  intensity  approaches  10  W/cm 
(Max  &  Estabrook  1979) .  A  large  fraction  of  the  absorbed  energy  is 
transferred  into  very  energetic  ("suprathermal" )  electrons  instead  of 
being  used  to  ionize  the  target.  Max  &  Estabrook  (ibid.)  postulate  that 


resonant  absorption  scales  as  IA  .  That  mechanism  is  dominant  above 
14  2 

10  W/cm  because,  as  the  charged  particles  become  more  energetic,  they 
exchange  a  relatively  smaller  fraction  of  their  kinetic  energy  on 
collisions  (Hora  1981) .  While  at  low  intensities  the  fraction  of 
incident  light  absorbed  by  the  target  is  high  (>80%)  at  both 
wavelengths,  at  intensities  of  order  1014  w/cm2  that  fraction  is 
typically  70%  at  0.35  urn  and  30%  at  1  ym  (Seka  et  al  1982).  IR  light  is 
substantially  less  effective  than  the  ratio  of  absorption  would  indicate 
because  the  suprathermal  electrons  serve  as  a  non-useful  channel  for 
deposition  of  the  incident  laser  energy. 

In  a  plasma  containing  mainly  Helium-like,  Hydrogen-like  and  bare 
ions,  the  limited  number  of  recombination  channels  is  an  advantage  for 
the  design  of  an  X-ray  amplifier  because  it  limits  the  waste  of  pump 
energy.  It  also  leads  to  simple  spectra  and,  with  proper  filtering,  the 
plasma  can  be  used  as  a  monochromatic  X-ray  source  (Conturie  &  Forsyth 
1980)  . 

4-1- a.  IR  illumination. 

Figure  4-1  shows  spectra  recorded  with  five  different  targets: 

Chlorine  (Z*17) ,  Calcium  (20) ,  Titanium  (22) ,  Iron  (26)  and  Nickel  (28) . 

15  2 

The  plasmas  were  produced  by  3x10  Watts/ cm  of  IR  light  (175  Joules  in 

700  psec  in  a  $  100  ym  focal  spot) .  The  pulse  duration  well  exceded  the 
ionization  equilibrium  time,  which  was  of  the  order  of  100  psec  (Whitney 
fc  Davis  1974) .  The  spectra  were  recorded  with  two  flat-crystal 
spectrographs  simultaneously,  at  45°  and  85°  from  the  laser  axis, 
equipped  with  collimating  slits  to  reduce  the  background  exposure  due  to 
suprathermal  electrons  generated  around  the  target.  The  most  prominent 


Figure  4-1  :  X-Ray  Line  Emission  from  Various 
Laser  Produced  Plasmas. 


group  of  lines  in  each  spectrum  includes  the  Helium-like  Is  -ls2p 
transition  (the  so-called  resonance  line,  labeled  w) ,  its  associated 
intercombination  lines  (labeled  x  and  y)  and  satellite  lines  arising 
from  dielectronic  recombination  to  doubly-excited  levels  in  the 
Lithium-like  ion  and,  for  Iron  and  Nickel,  in  the  Beryllium-like  ion  as 
well  (labeled  j  and  k) .  The  relative  intensity  of  the  resonance  line  in 
that  multicomponent  group  decreases  for  higher  atomic  numbers.  It 
becomes  harder  and  harder  to  ionize  the  atom  up  to  the  Helium-like  stage 
(not  to  mention  the  Hydrogen-like  stage)  and  the  approximately  constant 
electron  temperature  is  adequate  only  for  the  excitation  of  L-shell 
electrons  in  the  Lithium-like  and  Beryllium-like  configurations.  The 
overall  spectral  bandwidth  of  this  group  of  lines  is  of  the  order  of  1% 
or  less  so  that,  in  many  applications  (Frankel  &  Forsyth  1979),  the  fine 
structure  displayed  here  is  unimportant. 

The  nearest  spectral  component  to  this  group  is  the  ls-2p 
transition  (labeled  Lyman  a)  from  the  Hydrogen-like  ion  (and  associated 
Helium-like  satellites  lines) .  As  the  atomic  number  of  the  target 
increases,  the  intensity  of  L  drops  rapidly  compared  to  the  fall-off  in 
the  intensity  of  the  Helium-like  ion  resonance  radiation.  In  the  Iron 
spectrum  the  Hydrogen-like  emission  has  disappeared  altogether  and  the 
Helium-like  emission  is  sensibly  monochromatic  for  many  applications. 
Similar  spectra  were  recorded  on  targets  of  atomic  number  Z^14,  at  lower 
intensity  (see  Sec. 4-2) .  From  these  data,  and  using  the  film 
calibration  described  in  Appendix  A,  we  computed  the  number  of  X-ray 
photons  radiated  by  the  source.  Since  no  significant  difference  was 
found  between  the  two  directions  of  observation,  i.e.  45°  and  85°  from 


the  laser  axis,  we  derived  those  numbers  assuming  that  the  emission  from 
the  source  was  isotropic.  Figure  4-2  (lower  curve)  gives  the  fraction 
of  incident  laser  energy  converted  into  X-ray  energy  radiated  in  4tt  on 

the  Helium-like  resonance  line  group.  The  conversion  efficiency  is 

“3  ✓  -4 

about  10  for  Z*17  and  decreases  to  10  for  Z=28. 

For  elements  of  atomic  number  Z>15,  it  is  possible  to  isolate  that 
group  of  lines  by  taking  advantage  of  the  fact  that  its  wavelength  is 
slightly  longer  than  that  of  the  K-absorption  edge  of  the  solid 
material.  Therefore  a  thin  foil  of  the  target  material  may  be  used  to 
remove  the  other  line  components  present  in  the  spectrum.  An  example  is 
shown  in  Figure  4-3.  Under  identical  conditions  the  spectrum  of 
Chlorine  was  recorded  (a)  unfiltered  (same  spectrum  as  in  Fig.4-l-a) , 
(b)  through  0.5  mil  of  Saran,  (c2H2C*2*x'  ^  throu9h  1  mil  of  Saran. 

The  transmission  curve  of  Saran  is  given  in  Fig.3-3-a.  One  may  also 
notice  that  the  continuum  radiation  from  the  source  is  too  weak  to  be 
detected,  otherwise  one  would  see  a  difference  in  baseline  level  from 
one  side  of  the  Helium-like  resonance  line  to  the  other,  since  the 
transmission  through  the  Saran  changes  by  more  than  an  order  of 
magnitude. 

4-1-b.  UV  illumination. 

A  similar  series  of  experiments  was  conducted  using  UV  light  (by 
ourselves  in  our  experimental  vacuum  chamber,  the  "a- tank",  except  for 
Silicon  and  Nickel;  by  others  in  another  chamber,  the  "6-tank",  for 
those  two  elements  and  some  duplicates,  see  Table  A-l) .  Some  of  the 
spectra  are  shown  in  Figures  4-10  and  4-11.  Using  the  same  data 
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Figure  4-3  :  Monochromatization  of  Laser  Plasma 
X-Ray  Line  Emission  by  means  of 
Foil  Filtration. 


redaction  procedure  as  in  the  IR  case,  we  computed  the  X-ray  yield  on 
the  Helium-like  resonance  line  group  for  each  target  material.  As  is 
seen  in  Fig. 4-2  (upper  curve),  UV  irradiation  is  more  efficient  in 
producing  X-rays,  by  a  factor  (*10)  which  decreases  with  increasing  Z. 
This  observation  at  high  Z  is  believed  to  be  due  to  the  significant 
presence,  in  IR  laser  irradiation,  of  suprathermal  electrons  which  can 
contribute  to  the  excitation  of  high  energy  X-ray  lines  (Yaakobi  et  al 
1981-a).  The  main  reason  for  the  decrease  in  line  intensity  with 
increasing  Z  is  the  temperature  profile  in  the  target  (see  Fig. 1-1). 
Lines  of  higher-Z  ions  are  emitted  mostly  around  peak  temperature  near 
the  critical  layer  (where  the  laser  deposits  most  of  its  energy) , 
whereas  lines  of  lower-Z  ions  come  from  inside  the  critical  layer  where 
the  temperature  is  lower  and  the  density  higher. 

The  conclusions  drawn  from  these  experiments  are  the  following, 
a/  Laser-produced  plasmas  on  targets  of  atomic  number  z>15  can  be  used 
as  very  bright,  sub-nanosecond  point-sources  of  monochromatic  X-rays. 
Their  applications  as  a  probe  of  transient  phenomena  include:  (i) 
"Back-lighting"  of  micro-balloons  in  laser  fusion  (Key  et  al  1978).  An 
X-ray  source  synchronized  with  the  laser  driving  the  target  is  used  as  a 
probe  to  study  the  highly-compressed  imploded  fuel;  (ii)  X-ray 
diffraction  of  biological  samples.  Progress  has  been  made  in  that 
direction  with  the  recent  recording  of  the  patterns  from  crystal 
proteins  (Frankel  &  Forsyth  1982);  (iii)  Photopumping  leading  to  lasing 
in  the  soft  X-ray  region  (Jaegle,  et  al,  1981);  our  study  of  an  XUV 
amplifier  was  not,  however,  focused  in  that  direction;  (iv)  K-edge 


fine-structure  absorption  studies  (Yaakobi,  et  al,  1980). 


b/  Short  wavelength  lasers  couple  their  energy  to  the  target  more 
efficiently,  resulting  in  a  higher  X-ray  yield.  Using  UV  illumination 
instead  of  IR  leads  to  the  formation  of  a  higher  density  plasma,  at 
approximately  the  same  temperature,  a  significant  advantage  for  the 
design  of  an  X-ray  amplifier  based  on  three-body  recombination  (see 
Chapter  2).  In  such  an  amplifier  one  wants  to  form  a  large  reservoir  of 
bare  nuclei  (for  recombination  into  Hydrogen-like  ions) ,  or  of  H-like 
ions  (for  recombination  into  Helium-like  ions) .  With  the  lasers 
available  to  us,  this  sets  a  limit  for  the  atomic  number  of  the  target, 
Z<15. 


4-2.  Target  design  to  achieve  a  population  inversion. 

In  a  dense,  almost  fully  ionized  plasma  such  as  can  be  created  at 

the  surface  of  an  aluminum  target  (Z=13) ,  the  recombination  cascade  is 

dominated  by  collisions,  which  tend  to  populate  high-lying  quantum 

states  of  the  recombining  ion  (see  Chapter  2) .  In  particular  the 

three-body  recombination  rate  is  proportional  to  N2  while  the  radiative 

e 

recombination  rate  goes  as  Ng.  Also  the  probability  of  a  three-body 

event  is  strongly  dependent  on  the  electron  temperature  T  :  it  increases 

e 

more  rapidly  than  that  of  a  radiative  decay  as  Tg  cools  down  (Griem 

1964) .  Therefore,  the  idea  is  to  produce  a  high  density,  low 
temperature  plasma  from  an  initially  hot  plasma  (Gudzenko  6  Shelepin 

1965) .  One  may  rely  on  free  expansion  for  the  cooling  (Irons  &  Peacock 
1974,  Jacoby  et  al  1981)  but  it  is  accompanied  by  a  large  drop  in 
density,  greatly  reducing  the  possibility  of  significant  gain.  An 


alternate  solution  is  to  bring  the  plasma  in  contact  with  a  heat  sink 
(Bhagavatula  &  Yaakobi  1978):  a  foil  placed  a  few  hundred  microns  in 
front  of  the  target  allows  a  rapid  cooling  of  the  expanding  plasma  at  an 
early  stage  of  its  recombination,  by  heat  conductivity. 

We  pursued  the  latter  approach  and  devised  the  target  shown  on 
Figure  4-4-b.  The  foil  is  designed  as  a  slit  which,  if  properly 
dimensioned,  has  several  advantages:  (i)  large  surface  of  contact  with 
the  plasma,  without  vignetting  the  incident  laser  beam,  (ii) 
compatibility  with  a  point-focus  or  a  line-focus  illumination,  with  a 
direction  of  symmetry  ("x"  on  Fig. 4-4-b)  for  a  potential  X-ray 
amplifier,  (iii)  good  view  of  the  interaction  region  for  the  diagnostic 
instruments.  It  also  allows  one,  by  removing  one-half  of  the  slit,  to 
study  the  plasma-foil  interaction  separately  from  possible  nozzle 
effects  through  the  aperture. 

4-2-a.  Parametric  study  with  IR  illumination. 

At  1.054  ym,  the  target  (Fig. 4-4-b)  was  studied  over  a  parameter 
space  which  included  (i)  dimensions  of  the  foil  (gap  width  and  distance 
to  the  target,  (ii)  foil  and  target  materials,  (iii)  laser 
characteristics  (energy,  pulse  width) .  To  provide  a  baseline  for  the 
discussion  of  these  parameters,  we  present  on  Figure  4-5  a  spectrum 
(shot  #972)  which  illustrates  features  observed  in  much  of  the  data. 

The  laser  beam  (56  J,  600  psec,  f/3.5)  was  focused  on  an  aluminum 


slab.  Two  80  yn  foils  of  lead  (Pb,  Z*82)  forming  a  300  ym  slit  were 
placed  200  ym  in  front  of  the  slab.  The  spectrum,  recorded  with  a 
crystal  spectrograph,  is  a  direct  illustration  of  the  geometry  shown  in 
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igure  A-4  :  Scheme  for  an  X-Ray  Amplifier  /  Target  Geometry, 


Reference  Target  /  IR  Illumination 


Fiq.3-4.  The  set-up  was  such  that  the  magnifications  were  y'/y  =  2  and 

z'/z  =  3.3  (Al+^,  resonance  line)  to  4.2  (Al+*^,  Lyman  6).  The  scale 

on  Fig. 4-5  (and  all  through  Fig. 4-16)  refers  to  the  source  space. 

Figure  4-5-a  compares  a  densitometer  scan  of  the  dense  region  of 

the  plasma  (III)  with  those  of  the  two  spatial  integrations  (I)  and  (V) . 

Note  in  (V)  the  absence  of  satellite  lines  and,  compared  to  (III),  the 

smoother  decrease  along  the  H-like  Lyman  series  (L^  has  about  the  same 

intensity  in  the  two  spectra):  as  the  electron  temperature  T&  drops, 

recombination  into  the  higher  excited  states  is  favored.  In  (V)  the 
2  2 

Is  -ls3p.  Is  -ls4p  He-like  lines  and  even  the  L^ ,  L^  lines  present  ratio 
anomalies  large  enough  to  suggest  population  inversion  on  average  over 
the  expanding  plasma. 

The  last  discernible  lines  in  the  H-like  and  He-like  series  (ls-7p 
2 

and  Is  -ls6p  respectively)  give  an  indication  of  the  electron  density  N£ 
in  the  focal  region.  Using  the  Inglis-Teller  relation  (1939) ,  although 

it  may  not  be  valid  at  temperatures  higher  than  500  eV,  leads  to 

21  -3  21  -3 

N  =0.75X10  cm  in  one  case,  N  =1.7x10  cm  in  the  other.  This 
e  e 

21  -3 

corresponds  approximately  to  the  critical  density  (N  )  =10  cm  of  the 

e  c 
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plasma  at  X  =1.054  Um.  The  A1  (Is  -ls2p)  resonance  line  is  not  well 

ii 

resolved  from  its  associated  intercombination  and  satellite  lines  but  it 

+12 

appears  to  have  an  intensity  comparable  to  that  of  A1  (La) .  Thus  the 

electron  temperature,  on  the  basis  of  the  collisional-radiative  model 

which  applies  at  high  density  (Bates,  et  al,  1962-a) ,  is  estimated  to  be 

T  =750  eV. 
e 

The  diameter  of  the  plasma  in  the  focal  region  is  about  150  ym.  It 
may  be  larger  than  that  of  the  focal  spot  from  GDL  although  it  is  in 


aqceement  with  the  knife-edge  test  done  on  the  alignment  laser,  a  Nd : YAG 

which  goes  through  the  same  optics  as  GDL.  Thus  the  intensity  on  target 
14  2 

is  ^5x10  W/cm  .  Part  of  the  plasma  is  fully  ionized:  the  emission  on 

the  Lyman  series  when  the  expanding  plasma  reaches  the  foil  is  produced 

by  the  recombination  of  Al+^.  The  slight  asymmetry  in  the  wings  is  due 

to  a  miscentering  of  the  GDL  beam  through  the  foil.  Its  shape  is  the 

result  of  two  effects:  (i)  the  plasma  expands  around  each  foil  i.e. 

towards  the  center,  (ii)  the  plasma  emission  reflects  the  distribution 

of  the  electronic  heat  sink,  dense  near  the  foils  and  also  at  some 

distance  along  the  axis  where  its  components  merge. 

Figure  4-5-b  presents  scans  at  several  locations  of  the  shadowgram. 

Scan  #3,  which  represents  an  integration  over  the  region  y>400  ym,  shows 

+11  2  2 

an  intensity  inversion  between  the  Al  lines  Is  -ls3p.  Is  -ls4p  and 

+12 

possibly  between  the  Al  lines  Lg ,  (the  ambiguity  due  to  the 

Helium-like  continuum  precludes  a  definite  conclusion).  The  data 

reductions  from  the  spatially  resolved  spectrum  and  the  shadowgram  yield 

consistent  results:  away  from  the  foil  the  emission  on  the  Al+11  lines 
2  2  2 

Is  -ls3p.  Is  -ls4p.  Is  -ls5p  follows  an  exponential  decay  exp(-z/z  ) 

o 

with  zo=sl70  ym,  340  ym  and  430  ym  respectively.  Intensity  inversion  on 
Is3p/ls4p  and  Is3p/ls5p  occurs  200  ym  and  300  ym  away  from  the  foil 

while  it  is  never  observed  on  Is4p/ls5p.  Data  are  not  presented  for  L„ 

p 

and  Ly,  due  to  the  ambiguity  mentioned  earlier. 

4-2-a-l.  Dimensions  of  the  foil. 

Changes  are  observed  when  the  geometry  is  modified,  (i)  Localized 


cooling  with  very  little  deviation  of  the  plasma  expansion  is  observed 
on  shot  1547  (see  shadowgram.  Figure  4-6-a) ,  where  the  foil  is  a 


Figure  4-6  :  Lateral  Translation  of  the  Foil. 

Aluminum  (slab)  +  Magnesium  (half-plane  foil) , 

3.)  Shot  #547  (44  J,  600  psec,  IR)  ,  foil  at  the  edge  of  the  beam 
b)  Shot  #697  (63  J,  600  psec,  IR) ,  foil  moved  in. 


half-plane  of  magnesium  (thickness  %100  ym,  distance  to  the  slab  %  150 
ym)  positioned  just  out  of  the  incident  beam  (44  J,  600  psec) .  When  the 
foil  is  moved  in  (Fig.4-6-b,  shot  #697,  63  J,  600  psec),  the  plasma 
emission  follows  the  contour  of  the  foil,  (ii)  An  efficient  cooling  is 
achieved  when  the  foil  is  200  ym  to  300  ym  wide,  about  200  ym  from  the 
slab.  The  slit  is  a  compromise  between  obstruction  of  the  laser  beam 
and  proximity  to  the  plasma  expansion. 

4-2-a-2.  Foil  and  target  materials. 

The  choice  of  foil  material  does  not  strongly  affect  the  result. 

One  shot  (#978,  54  J,  600  psec,  see  Figure  4-7-b)  was  taken  with  a  slit 
made  of  one  half-plane  of  magnesium,  the  other  of  lead  for  direct 
comparison.  No  clear  spatial  asymmetry  is  apparent,  which  is  consistent 
with  what  is  observed  on  shots  taken  with  either  type  of  foil  alone.  In 
many  subsequent  experiments,  we  used  a  foil  of  Titanium  (z*22)  because 
(i)  it  does  not  lead  to  spectral  overlap  with  the  lines  from  the 
original  plasma  (Z<15) ,  (ii)  it  has  a  relatively  high  melting  point, 
1668°C,  and  does  not  produce  a  blow-off  damaging  to  the  optics  (as  a  Pb 
foil  does) . 

Other  target  materials  were  used:  Magnesium  (Z«12)  and  Silicon 

(Z»14) .  The  results  are  shown  in  Fig.4-7-a  (shot  #4123,  54  J,  600  psec) 

* 

and  4-7-c  (shot  #1740  ,  20  J,  600  psec),  respectively.  As  before,  there 
is  reemission  of  the  original  plasma  as  it  reaches  the  vicinity  of  the 
foil,  and  population  inversion  occurs  in  the  downstream  expansion. 

That  set  of  data  allows  us  to  put  aside  models  based  on  resonance 


mechanisms,  such  as  photo-pumping  (Bhagavatula  1980)  and  charge  transfer 
(Elton  et  al  1979).  Both  are  strongly  material  dependent.  In 
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Figure  4-7  :  Targets  of  Different  Materials  (Slab  +  Foil). 

a)  Shot  #4123  (54  J,  600  psec,  IR) ,  Mg  +  Ti 

b)  Shot  #  978  (54  J,  600  psec,  IR)  ,  Al  +  Mg/Pb 

c)  Shot  #1740* (20  J,  600  psec,  IR) ,  Si  +  Mg 
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photo- pumping  there  must  be  resonant  absorption  of  the  pump  radiation,  a 
condition  not  satisfied  in  most  combinations  of  target  and  foil 
materials  mentioned  above.  In  charge  transfer  it  was  suggested  (ibid.) 
that  Neon-like  ions,  which  have  a  closed  shell  and  are  expected  to  be 
long-lived  even  in  a  transient  plasma,  could  transfer  an  electron 
selectively  to  an  excited  level  of  hydrogenic  or  Helium-like  ions  of  the 
same  species.  This  would  require  a  match  between  the  binding  energies 
of  the  electron  in  the  two  configurations,  with  some  energy  defect  to 
insure  an  exothermic  exchange.  These  conditions  could  be  met  in  Al  or 
Mg  but  not  in  Si. 

To  confirm  the  hypothesis  of  a  lowering  of  the  plasma  temperature 

near  the  foil,  leading  to  enhanced  recombination,  we  used  a  target  of 

alumina  (Al^^)  instead  of  aluminum .  Figure  4-8-b  shows  the  oxygen 

* 

(Z«8)  spectrum  recorded  on  shot  #1108  (30  J,  600  psec.  Mg  foil).  In 

the  focal  region  (which  is  within  the  field  of  view  of  the  spectrograph, 

+12 

as  indicated  by  the  recording  of  the  Al  Lyman  ot  line  in  third  order) , 
the  temperature  is  too  high  for  fully-stripped  oxygen  ions  to  recombine. 
Most  of  the  line  emission  occurs  near  the  foil.  Without  foil  the 
emission  is  more  diffuse,  because  it  occurs  during  the  plasma  expansion, 

and  only  the  space- integrated  spectrum  can  be  recorded  (see  Fig.4-8-a, 

* 

shot  #970  ,  20  J,  600  psec). 

4-2-a-3.  Laser  characteristics. 


Most  of  the  data  presented  so  far  were  recorded  on  IR  shots  of 
energy  ^50  J  and  duration  ^>600  psec.  While  plasma  reemission  near  the 
foil  is  observed  at  laser  energies  as  low  as  6  J,  the  observation  of  a 
population  inversion  in  the  expanding  plasma  of,  say,  Al+^  requires  at 


vV> 


least  20  J.  Below  20  J,  the  pump  is  not  strong  enough  to  create  a  large 
reservoir  of  bare  nuclei  or  hydrogenic  ions.  Above  50  J,  most  of  the 
energy  increment  contributes  to  the  production  of  suprathermal 
electrons,  which  do  not  participate  in  the  formation  of  that  reservoir. 

A  few  spectra  were  recorded  with  shorter  pulses  (~150  psec) .  The 
familiar  pattern  of  X-ray  emission  in  the  focal  region  and  near  the  foil 
was  observed  again  but  without  population  inversion.  The  onset  of 
suprathermal  electrons,  which  depends  on  the  laser  intensity  on  target, 
occurs  at  a  lower  energy  level  than  in  long  pulse  operation.  The  pump 
operates  for  a  shorter  time,  without  improved  efficiency  or  any  other 
advantage. 

4-2-a-4 .  Other  targets. 

Two  variations  of  the  step  target  were  studied,  with  the  same  idea 
of  enhanced  recombination  in  mind.  Both  were  designed  (i)  to  keep  the 
plasma  density  high  by  containing  its  expansion,  (ii)  to  improve  the 
cooling  by  increasing  its  area  of  contact  with  the  heat  sink. 

The  first  was  a  channel  (Reintjes  et  al  1978,  Silfvast  et  al  1980)  , 

i.e.  a  thick  foil  whose  distance  to  the  target  was  zero.  Figure  4-9-a 
* 

(shot  #1826  ,  51  J,  600  psec)  shows  the  spectrum  obtained  with  an  Al 
target  and  a  Mg  channel  (300  ym  thick,  250  ym  wide).  The  second  was  a 


tilted  target  (see  Fig.4-4-c).  Even  though  the  expansion  is  expected  to 
be  mainly  in  the  direction  of  the  incident  laser  (Hora  1981) ,  the  design 
lengthens  the  interaction  time  between  the  plasma  and  the  foil. 

Fig.4-9-b  (shot  #3091,  66  J,  750  psec)  shows  the  spectrum  obtained  with 
an  Al  target  and  a  Mg  foil. 
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Figure  4-9  :  Alternate  Target  Designs, 

Aluminum  (slab)  +  Magnesium  (foil), 

a)  Shot  #1826* (51  J,  600  psec,  IR)  ,  Channel  Target 

b)  Shot  #3091  (66  J,  750  psec,  TR) ,  Tilted  Target. 
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In  both  cases  a  population  inversion  is  observed  between  the  levels 
n*4  and  n«3  of  the  Al+^  ion.  The  inversion  is  of  a  magnitude  and  at  a 
location  in  the  expansion  similar  to  that  observed  in  shot  #972 
(Fig. 4-5).  Hence  it  appears  that  what  counts  to  produce  enhanced 
recombination  is  the  presence,  at  an  early  stage  of  the  expansion,  of  a 
foil.  The  geometry  of  the  foil  determines  the  path  of  the  plasma  flow, 
but  the  recombination  is  not  influenced  by  it  being  a  barrier  (Fig. 4-6), 
a  nozzle  (Fig. 4-5)  or  a  channel. 

4-2-b.  UV  illumination. 

Experiments  on  step  targets,  similar  to  those  described  in 
Sec.4-2-a,  were  done  with  the  frequency- tripled  GDL  laser  (0.351  ym)  to 
take  advantage  of  the  superior  coupling  of  short  wavelength  light  to  the 
target.  Figure  4-10  shows  the  spectrum  recorded  when  the  UV  laser  (shot 
#2789,  20  J,  500  psec)  is  focused  in  a  4>  100  ym  focal  spot  at  the 
surface  if  an  aluminum  slab.  Two  100  ym  foils  of  titanium,  forming  a 
200  ym  slit,  were  placed  200  ym  in  front  of  the  slab.  As  before, 
aluminum  X-ray  emission  is  observed  mainly  in  the  focal  region  at  the 
surface  of  the  target  and  along  the  contour  of  the  foil,  with  a  more 
diffuse  tail  as  the  plasma  expands  beyond  the  foil.  As  is  apparent  on 

the  densitometer  trace,  there  is  an  intensity  inversion  at  distances 

2  2 

600  ym  into  the  expansion  between  the  Is  -ls3p  and  Is  -ls4p  lines  in 
Al+^,  showing  the  possibility  of  gain  at  129.7  A  on  33D  -  43F.  Without 
foil,  emission  strong  enough  to  be  recorded  occurs  only  in  the  focal 
region  where  high  density  and  temperature  adequate  for  line  excitation 
can  be  found  along  the  gradient  from  the  hot  corona  to  the  cold  target 
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mass.  The  expanding  plasma,  initially  at  high  temperature,  is  too 

dilute  for  measurement  by  the  time  it  has  cooled  down  to  the  point  where 

the  rate  of  ion  recombination  (and  line  radiation)  is  largest. 

An  attempt  was  made  to  scale  the  experiment  up  to  higher  atomic 

numbers.  The  UV  laser  was  focused  on  a  target  of  similar  geometry,  but 

whose  slab  and  foil  were  of  Titanium  (Zs22)  and  Magnesium,  respectively. 

Once  again  the  choice  of  foil  material  was  mainly  based  on  avoiding 

spectral  overlap.  Figures  4-11-a  and  4-11-b  show  the  spectra  recorded 

* 

when  the  UV  laser  is  focused  on  such  a  target  (shot  #3712  ,  37  J,  1 

* 

nsec),  and  on  a  flat  Titanium  target  without  foil  (shot  #3713  ,  32  J,  1 

nsec),  respectively.  Because  of  the  higher  atomic  number  compared  to 

Al,  it  is  harder  to  achieve  the  corresponding  ionization  state.  Thus, 

the  smaller  reservoir  of  hydrogenic  and  fully-stripped  ions  leads  to  a 

much  reduced  level  of  recombination  into  the  helium-like  stage  when  the 

plasma  reaches  the  foil.  Line  reemission  is  still  apparent  on  the 

2  i  3 

Helium-like  resonance  transition  group  (Is  -ls2p  "p  6  P)  and,  very 
2  2 

weakly,  on  the  Is  -ls3p  and  Is  -ls4p  transitions.  There  is  evidence  of 

enhanced  recombination,  as  can  be  seen  from  the  significant  increase  of 
2  2 

the  ratio  Is  -ls4p/ls  -ls3p  in  the  expansion,  to  ^1  about  500  ym  from 
the  target  surface.  The  result,  while  marginal,  illustrates  the 
possibility  of  scaling  towards  shorter  wavelengths  when  larger  pump 
lasers  are  available.  The  transition  Ti+^® (3^D-4^F)  is  at  43  A. 

4-3.  Design  of  an  ASE  device. 

An  XUV  amplifier  requires  an  elongated  medium  (i)  to  maximize 
the  length  of  amplification,  (ii)  to  permit  an  unambiguous  diagnostic  of 


gain  from  the  spatial  distribution  of  the  amplified  spontaneous  emission 
(Slaymaker  1978) .  That  requirement,  in  schemes  involving  traveling 
waves  (Chapline  &  Wood  1974,  Duguay  1976),  is  an  essential  part  of  the 
design.  In  a  steady-state  type  of  pumping  such  as  ours,  the  best 
approach  is  to  use  cylindrical  optics  to  form  a  line-shaped  plasma.  The 
optics  are  described  in  Appendix  C.  However,  as  the  focal  area  gets 
larger  with  a  cylindrical  lens,  the  intensity  on  target  decreases  and 
one  needs  to  use  materials  of  lower  atomic  number  (Z)  to  form  a  fully 
ionized  plasma.  In  addition,  as  Z  decreases,  the  recorded  spectrum  is 
simpler  and  the  detection  of  the  relevant  spectral  lines  is  made  easier. 

4-3-a.  Aluminum  (Z=13)  and  Magnesium  (12) . 

Experiments  with  IR  light  were  done  using  a  double  aspheric  lens 
which  produced  a  1  mm  line-focus  of  aspect  ratio  '''10:1.  Figure  4-12-a 
shows  a  pin-hole  picture  of  the  X-ray  emission  (5  k<\<9  %)  from  the 
plasma  formed  at  the  surface  of  an  aluminum  target  (shot  #3311,  78  J, 

600  psec,  Mg  foil).  Figure  4-12-b  shows  a  simulation  with  the 
ray-tracing  code  LINFOC  (see  Appen.C) .  However  the  data  taken  with 
crystal  and  grating  spectrographs  on  that  shot  (as  well  as  others)  point 
to  several  difficulties: 

14  2 

a)  the  intensity  on  target  (1=10  W/cm  )  is  too  low  to  create  a 
population  inversion  in  hydrogenic  or  helium-like  aluminum,  at  least  not 
one  which  could  be  recorded  on  a  time-integrated  basis. 

b)  target  oxidation  prior  to  the  onset  of  vacuum  leads  to  the 
recording,  in  the  100  &  region,  of  an  oxygen  spectrum  superimposed  on 


the  spectrum  of  interest. 
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Figure  4-12  : 


-  Pin-hole  picture. 

Aluminum  (slab)  +  Magnesium  (foil) 

-  Notice  the  X-ray  emission  near  the  foil 
on  the  left  side  of  the  slit. 

-  Vertical  dimension  reduced  by  fl 
(observation  at  45°) . 
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a)  Shot  #3311  (78  J,  600  psec,  IR) 

b)  LINFOC  simulation. 


c)  laser-transition  candidates,  such  as  Al+*3 (33D-43F)  at  129.7  A, 
cannot  be  resolved  because  they  overlap  with  strong  lines  from  other 
ionic  stages  of  A1  or  0  (see  Chap. 3). 

All  these  problems  can  be  solved  by  using  UV  light  and  a  lower-Z 
target  (see  next  section).  It  may,  however,  be  possible  to  build  a 
device  and  do  the  diagnostic  without  significantly  changing  Z.  Problem 
b  can  be  solved  by  having  a  pre-shot  on  the  target  (while  protecting  the 
foil  and  closing  the  aperture  of  the  spectrographs) ,  to  remove  the  thin 
layer  of  surface  oxide.  An  illustration  of  that  idea  is  shown  on  Figure 
4-13,  which  is  a  comparison  of  the  spectrum  from  a  flat  aluminum  target 
with  and  without  a  "clean-up"  shot  beforehand.  The  oxygen  lines 
disappear  almost  completely  when  that  procedure  is  used  (shot  #4021) . 
Problem  c  can  be  solved  by  using  a  magnesium  target  (Z=12) .  Figure  4-14 
shows  that  the  line  Mg+1° (33D-43F)  at  154.4  &  can  be  resolved  from  its 
neighbors.  The  full  spectrum  in  the  range  45  A-200  A  was  recorded  on 
film  (see  Chap. 3),  the  wavelength  scale  was  determined  with  the  code 
GRASPEC  (see  Appen.B) ,  and  the  lines  were  found  in  good  agreement  with 
listed  wavelengths  (Kelly  &  Palumbo  1973) .  Most  of  the  features  around 
154  X  are  actually  second  order  lines.  The  precise  correspondence 
between  the  first  and  second  order  spectra  give  added  confidence  in  the 
identification  of  the  transition. 

4-3-b.  Fluorine  (Z«9) . 

When  the  focal  spot  of  the  laser  is  extended  from  a  point  to  a 


line,  the  intensity  on  target  with  currently  available  systems  is  such 
that  only  lower-Z  materials  are  suitable  for  a  population  inversion 


Removal  of  the  Oxydized  Surface  Layer 
of  an  Aluminum  Target. 
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RECORDING  OF  3d  -  4f  TRANSITION  (154.4&) 
IN  HELIUM-LIKE  MAGNESIUM 


3  *  (30  J,  600  psec) 
1.054  pi  m  Light 
Flat  Target 


Mg  VIII  j-—  2p3-  2p2(3P13d  2a2p2  -  2s2p  [ 1  3d 


PP]  l3P] 


First  Order 
ol  Diffraction 


ZT9 

*  . 

k-KT 

> 

'  -  «~ 

> 

rv 

.  >  V 

-  6 

k  0  t 

V 

* 

e 

*•*  >■ 

a  ; 

.  r 

*  • 

-  •  * 

_ at 

"V 


Second  Order 


VI 


154.4  A 

_J _ i 

154 


v  y 


Figure  4-14 


study.  Figure  4-15-a  shows  the  fluorine  spectrum  recorded  when  the  UV 
* 

laser  (shot  #4471  ,  15  J,  500  psec)  is  focused  on  a  100  pm  x  2  m  area 

of  teflon  ((CF2^n^*  Two  t*tanium  foils,  100  ym  thick  and  forming  a  200 

ym  wide  slit,  were  placed  200  ym  in  front  of  the  target.  For 

comparison,  Fig.4-15-b  shows  the  spectrum  recorded  when  the  foils  are 

* 

removed  (shot  #4475  ,  20  J,  500  psec). 

* 

In  shot  #4471  ,  the  plasma  expansion  is  very  directional,  planar 

instead  of  diverging.  As  is  apparent  on  the  film  and  densitometer 

traces,  the  transition  Lyman  g  (ls-3p)  is  much  more  intense  than  Lyman  a 

(ls-2p)  almost  everywhere,  on  a  time-integrated  basis.  There  are 

several  caveats:  (i)  the  heat  sink  partially  blocks  the  line  of  sight 

of  the  spectrograph  for  the  radiation  emitted  near  the  target,  for  X>14 

A,  (ii)  the  entrance  window  of  the  spectrograph  (0.75  ym  aluminum) 

transmits  about  15%  more  at  12.64  A  (Lg)  than  at  14.98  X  (La) ,  (iii) 

the  film,  RAR-2491,  is  more  sensitive  at  the  shorter  wavelength 

(Benjamin  et  al  1977).  However,  even  after  taking  those  factors  into 

4*8 

account,  there  exists  a  strong  intensity  inversion  between  the  two  F 
lines  through  most  of  the  expansion,  showing  the  possibility  of  gain  at 
80.91  X.  There  is  also  some  evidence  of  an  inversion  between  Ly  (ls-4p) 
and  Lg  late  in  the  expansion  (  >1  mm  from  the  target  surface) .  Gain 
would  be  at  231.3  X  on  3d-4f  and  at  59.95  X  on  2p-4d. 

An  attempt  was  made  to  record  these  transitions  directly  with  a 
grating  spectrograph.  Only  a  direct  observation  will  lead  eventually  to 


at  teflon  target,  spatially  integrated  spectrum 
#4518* ( 30  J,  1  nsec,  UV,  line-focus), 
tanium  foil  added,  spatially  resolved  spectrum. 
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a  conclusive  demonstration  of  XUV  amplification.  If  amplification 
occurs,  the  radiation  emerging  from  the  ends  of  the  elongated  plasma 
should  be  more  intense  than  that  emitted  in  other  azimuths,  because  the 
photons  travel  a  longer  distance  through  the  medium.  Therefore,  the 
angular  distribution  of  the  amplified  spontaneous  emission  (ASE)  should 
be  peaked  along  or  near  the  end  direction  (Slaymaker  1978) .  A 
meaningful  diagnostic  requires  good  discrimination  between  the  angularly 
dependent  signal  and  other  types  of  radiation.  The  transition  of 
interest  must  be  recorded  (i)  on  a  single-shot  basis,  to  avoid  problems 
of  reproducibility  in  illumination  conditions  and  target  design,  (ii) 
with  spatial  resolution  along  the  plasma  expansion,  to  isolate  the  "gain 
region",  (iii)  free  of  overlap  with  other  spectral  components  or  strong 
continuum.  Time  resolution  may  also  be  required,  due  to  the  transient 

nature  of  the  gain  (see  Chap. 5).  Figure  4-16  shows  the  spectrum 

* 

recorded  on  shot  #4518  (30  J,  1  nsec,  same  target  and  focusing 

* 

conditions  as  shot  #4471  ),  with  an  80  ym  spatially  resolving  slit.  A 

* 

reference  spectrum  (shot  #2906  ,  10  J,  600  psec,  point  focus  if  100  ym), 
taken  without  spatial  resolution,  is  also  presented.  The  spectral 
resolution  in  both  cases  is  A 0.2  &.  Hydrogenic  (F  IX)  and 
Helium-like  (F  XIII)  emissions  are  recorded  up  to  distances^  1  mm  and  ~2 
mm  from  the  target  surface,  respectively.  Transitions  from  lower  ionic 
species  are  still  observed  at  the  edge  of  the  instrumental  field  of 
view,  3  mm  from  the  target  surface.  The  Balmer  ot  line  (F  IX,  2-3)  is 
'well  isolated  from  other  features  and  stong  enough  to  be  a  suitable 
candidate  for  an  actual  study  of  the  ASE.  The  value  of  the  gain  x 
length  product  is  somewhat  uncertain  due  to  the  limitations  of  our 


model,  as  will  be  seen  in  the  next  chapter  (Sec.5-3-c).  Indications  are 
that  it  is  still  below  1,  even  on  an  instantaneous  basis,  and  that  an 
angular  measurement  of  the  ASE  would  not  be  conclusive.  Slaymaker 
(ibid.)  calculates  that,  depending  on  the  radiation  environment, 
i.e.  the  noise  level,  a  gain  x  length  product  >5  is  required  for 
detection.  It  is,  however,  a  very  promising  direction  of  research  and 
Fig. 4- 16  shows  that  the  diagnostic  capability  is  in  place,  ready  to  be 
used  when  higher  gain  is  achieved. 

In  this  chapter,  we  have  shown  that  (i)  short  wavelength  lasers 
couple  their  energy  to  the  target  more  efficiently,  (ii)  targets  of 
atomic  number  Z>15  can  be  used  as  sub- nanosecond  point  sources  of 
monochromatic  X-rays,  (iii)  targets  of  atomic  number  Z<15  can  be  used 
for  the  design  of  an  X-ray  amplifier  based  on  enhanced  recombination  of 
an  expanding  plasma,  with  good  scaling  prospects  when  larger  systems  are 
available,  (iv)  the  Balmer  o t  line  of  fluorine  (80.9  &)  is  a  prime 
candidate  for  direct  gain  measurement  in  the  not-too-distant  future. 


Computer  Modeling 

The  experiments  described  in  Chapter  4  were  modeled  with  the 
1-dimensional  laser  fusion  code  LILAC  (Goldman  1973,  Delettrez  1978-82). 
The  code  computed  the  hydrodynamic  motion  of  the  plasma  while  providing 
temperature  and  density  of  both  ions  and  electrons  at  each  time  step  and 
for  each  numerical  cell.  These  parameters  served  as  input  for  our 
subroutines  modeling  the  atomic  rate  equations,  radiation  transport  and 
plasma  cooling  relevant  to  the  present  study.  This  chapter  gives  some 
details  about  the  program  (file  name  CYLINE)  containing  the  subroutines. 
A  dummy  calling  program  (YRATE)  was  written  to  run  simulations  for 
predetermined  expansion  characteristics  (N^  T*)  and  cooling  times. 
Full-scale,  time-dependent  simulations  with  LILAC  (in  place  of  YRATE) 
are  compared  with  the  experimental  data. 

5-1.  Program  CYLINE. 

CYLINE  is  based  on  the  calculations  by  McWhirter  &  Hearn  (1963)  for 
the  recombination  of  a  fully-ionized  plasma  into  an  hydrogenic  plasma, 
as  reviewed  in  Chapter  2.  We  describe  here  some  of  the  specific 
features  of  our  program. 

5-1- a.  Rate  equations. 

The  treatment  of  McWhirter  &  Hearn  (ibid.)  was  also  used  for  the 


Helium-like  ions  with  (i)  an  apparent  nuclear  charge  Z-l  instead  of  Z 
for  the  reduced  electron  temperature  and  density,  (ii)  the  assumption  of 


a  uniform  population  density  among  states  of  same  principal  quantum 

number,  (iii)  a  proper  count  of  the  degeneracy  factor  for  each  energy 

level.  Outside  the  range  of  validity  of  the  calculations  (0.35 
2  2 

Z<Te(eV)<22  (Z-l)  ,  no  practical  boundaries  for  N^) ,  the  populations  of 
the  excited  states  were  not  computed.  It  is  not  a  serious  limitation 
because  (a)  in  the  upper  range  of  temperature  the  plasma  contains 
mainly  bare  nuclei,  (b)  in  the  lower  range  most  hydrogenic  and 
helium-like  ions  are  in  their  ground  state.  For  lower  stages  of 
ionization,  we  used  a  formula  quoted  by  Stratton  (1965) ,  which  applies 
to  a  plasma  in  coronal  equilibrium: 


log  (N .  /N . )  -  8.15  +  log  T  -31ogx.  -X-/2*3T  +log(£  /n  )  (5-1) 

x+1  i  e  i  l  e  oo 


"where  and  Ni+1  are  the  relative  numbers  of  ions  of  states  of 


ionization  i  and  i+1,  (...)  X^  is  the  energy  [in  eV]  required  to  ionize 


the  ith  ionization  state,  (...)  n  is  the  total  quantum  number  of  the 

o 


ground  state  of  the  ith  ionized  atom,  (...)  £  is  the  number  of 

n 

o 


electrons  in  the  n  th  shell,  (...)  and  T  is  the  electron  temperature 
o  e 


(in  eV]"  (ibid.).  Coronal  equilibrium,  i.e.  the  regime  of  collisional 
excitation  balanced  by  radiative  recombination,  applies  when  (ibid.) 


N  <  3  1013  X?  <kT  )1/2 
e  i  e 


(5-2) 


Since  the  rapid  laser  heating  produces  a  high-energy  plasma,  L-shell 
recombinations  become  significant  only  after  the  plasma  has  expanded  and 
cooled.  Condition  5-2  then  holds  true  and  justifies  the  use  of  Eq.5-1. 

It  is  important  to  adapt  the  rate  equations  to  the  constantly 
changing  temperature  and  density  of  the  plasma.  The  rates  of  three-body 
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and  radiative  recombination  and  ionization  (Elton  1970)  were  used  to 

determine  the  plasma  equilibrium  time,  x  .  x  was  set  equal  to  the 

eq  eq 

fastest  of  these  rates  and  was  compared  to  the  time  step,  X  of  the 

h 

hydrodynamic  calculation.  If  T^<Teq»  the  populations  (N^J^  calculated 
at  step  j  do  not  change  rapidly  enough  to  reach  their  new  equilibrium 
value  eg  at  step  j+1.  The  actual  populations  (N^)  j+1  were 

computed  as  follows: 


‘Vj.l  -  (Vj  +  '  '‘Vju’eq'  (»i)j)U-e*P(-Th/Teij)l  ,5-3) 

Whenever  necessary,  the  population  distribution  was  renormalized  to  the 
total  number  of  ions  in  each  numerical  cell. 


5-1-b.  Radiation  transport. 

CYLINE  calculates  the  number  of  photons  received  by  a  detector 
placed  at  90°  from  the  laser  axis  (see  Figure  5-1) .  The  "detector" , 
which  simulates  a  crystal  spectrograph,  provides  a  10  ym  spatial 
resolution  along  the  direction  of  expansion  (that  resolution  could  be 
changed  easily  if  necessary) .  In  a  dense  plasma,  strong  reabsorption  of 
line  radiation  may  occur.  It  affects  not  only  absolute  photon  counts 
but  also  line  ratios  because  it  is  proportional  to  the  absorption 
oscillator  strength  (f)  of  the  transition.  Let  us  assume  a  uniform 
plasma  of  length  1  and  absorption  coefficient  (see  Fig. 5-1-b).  The 
intensity  measured  at  one  end  of  the  plasma  is  (Elton  1970) 

\  m  11  -exP<“Kv1)]  (5-4) 

where  £v  i®  the  emission  flux  density  (in  the  absence  of  absorption 
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EQUIVALENT  TARGET  IRRADIATION  GEOMETRY 
USED  IN  ONE-DIMENSIONAL  (LILAC)  SIMULATION 


Figure  5-1 


(<v*0) ,  I  *  evl).  Hence  the  transmission  factor  for  the  radiation  at 
frequency  v  is 


t  -  tl-exp(-<vl)]/<vl  =  [l-exp(-Tv)]/Tv  (5-5-a) 

where  we  define  the  optical  depth,  1.  Eq. 5-5-a  has  to  be 

integrated  over  the  normalized  spectral  profile  of  the  transition 

»  3 

(  :  dv  S  =N,  density  of  absorbing  ions  per  cm  ): 

-oo  V 

too  00 

=  f  dv  S  {[  (1  -  exp  (-T  )  ]  /T  }  /  I  dv  S„  (5-5-b) 

V  —  oo  V  V  V  —oo  V 


Since  Tv  is  proportional  to  S  , 


OO  00 

I  dv  [1  -  exp  (-T  )  ]  /  /  dv  T 

•oo  V  —oo  V 


(5-5-c) 


Let  us  consider  the  two  basic  cases  of  Lorentzian  and  Doppler  spectral 
profiles. 


5-1-b-l.  Opacity  for  Lorentzian  lines.  If  vq  is  the  frequency  at 
the  center  of  the  line  (optical  depth  T),  and  Av  the  fwhm,  (Svelto  1982) 


T  =  T(Av/2)2/[(v-v  ) 2+ (Av/2) 2] 
V  o 


(5-6) 


In  the  numerator  on  the  right-hand  side  of  Eq. 5-5-c,  the  exponential  can 
be  written  as  a  series.  The  term  by  term  integration  leads  to  (5-7) 


/dv  II  -exp(-T  ) )  «  I(-l)n+1(Tn/n!) (Av/2)2n  /dv/[(v-v  )2+(Av/2)2]n 

•°°  v  n*l  —oo  O 


From  the  identity  (Gradshteyn  &  Ryzhik  1980,  p.292) 

1—  2n  n  .  ■  j.%  «  >  -  >  _  >  i  n  <  i 


■  a 


/  dx/  (x2  +  a2)n 


ra/2)  r(n-i/2)  /  r(n) 


(5-8) 


it  follows  that  the  right-hand  side  of  Eq.5-7  is  equal  to 


tt (Av/2)  {  ?.(-l)n+1  <r(n-i)  /  r<4>>  Tn/(n-l)!  n'} 
n=i 


(5- 


Pinally,  after  dividing  by  /dv x  -  tt(Av/2)t,  we  get 

-oo  V 


fT  v  «  1  +  ?,(r(n+i)  /  r<A))  (-T)n/  n!  (n+1) ! 

Lorentz  n=l  2  2 


(5' 


The  series  converges  because  (i)  consecutive  terms  are  of  opposite 
sign,  (ii)  the  ratio  of  consecutive  terms  is  (n-0.5)  fc/n(n+l) <  l  for 
large  values  of  n  (n>t). 


5-l-b-2.  Opacity  for  Doppler  lines.  In  this  case  (Svelto  1982) 


*  t  exp- [A(v-vo) ] ' 


(5 


where  A  is  a  constant  which  need  not  be  specified,  as  will  be  apparent 
in  the  final  result.  The  same  series  expansion  of  the  exponential  in 
Eq.5-5-c  leads  to 


/  dv  [1  -  exp(-T  )  ]  =  {  Z  (-l)n+1  xn/n.'/n}  7  dv  exp-[A(v-v  )]2  (5 

.oo  v  n=l  ° 


After  division  by  /dv  x  , 


t 


(-1/T)  Z  (-T)n/n!v^  =  Z  t 
n=l  n-1  n 


'Doppler 

The  series  converges  for  the  same  reasons  as  the  previous  one  did. 


(5 


The  program  CYLINE  assumes  Doppler-broadening  of  the  spectral 
lines.  Ion  temperatures  T^  in  the  range  0.1-1  keV  produce  linewidths 

.A 

(Av/v)  ,  of  order  1.5-5x10  in  aluminum  (see  Eq.3-1).  Stark 
Doppler 


’w 
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broadening,  which  would  produce  a  Lorentzian  line  profile,  occurs  at 
densities  higher  than  those  of  the  expanding  plasmas  studied  here 
(Yaakobi  et  al  1981-b) .  Also,  since  it  is  proportional  to  the 
difference  of  the  squared  total  quantum  numbers  of  the  upper  and  lower 
levels,  it  should  be  larger  for  transitions  from  high- lying  states. 
That  effect  is  not  observed  on  any  of  the  data  presented  in  Chapter  4. 

The  optical  depth  at  the  center  of  the  Doppler  line  profile  is 
given  by  (Griem  1964,  p.196): 


T 


1.1  10‘16  X  K-<1>  (M/T^) 1/2  f 


(5-14) 


where  M  is  the  ion  mass  (amu) ,  T^  the  ion  temperature  (eV) ,  X  the 

wavelength  (A)  and  < 1>  (cm)  the  average  length  of  the  plasma  in  the 

direction  of  observation.  Refering  to  Pig. 5-1,  (i)  in  a  point-focus 

geometry,  <  1>*  (8/3tt)  ($/r )  z  where  <j>  is  the  diameter  of  the  focal  spot,  r 

the  radius  of  the  sherical  target  used  in  the  simulation,  z  the  radial 

distance  and  8/3rr  an  averaging  factor  over  the  circular  cross-section  of 

the  expanding  plasma,  (ii)  in  a  line-focus  geometry,  < l>*Lf +0.5 (z-r) 

where  is  the  length  of  the  focal  spot  and  0.5  a  factor  chosen  to 

simulate  an  expansion  in  a  cone  of  half-angle  0=14°.  in  order  to 

compute  Eq.5-13,  we  established  (with  a  little  code,  file  name  OPACITY) 

«£  "1 

an  asymptotic  relation  for  large  values  of  T,  *^p0pp]_er~  (1*5+0.40 
For  t<30,  CYLINE  computes 


"l  t  +1/2  t  * 
n=l  n  n  +1 

* 

where  n  »3  Integer (T)  +  10.  For  x>30,  it  computes 


( 5-15-a) 


f  -  £  13. 5/(1. 5+0. 4  T) 


( 5-15-b) 


5-1-c.  Plasma  cooling. 

The  effect  of  the  foil  in  front  of  the  target,  which  acts  as  a  heat 

sink,  is  treated  by  a  simple,  spatially  dependent  energy  loss  term 

proportional  to  Tg ( 1-exp (-t/xc) )  where  t  is  the  instantaneous  time  step 

of  the  calling  program  (LILAC  or  YRATE) ,  and  t  is  a  "characteristic 

c 

cooling  time"  used  as  free  parameter.  In  CYLINE  itself,  that  loss  is 
written  in  the  form  of  a  temperature  drop: 

-  AT  /T  =  1  -  exp{-(t/T  )  exp- r  (2z-z  -z  )/(z  -z  )  ]4}  (5-16) 

e  e  c  2121 

where  z^  and  z^  are  the  radial  coordinates  of  the  foil  boundaries.  The 
spatial  modulation  was  introduced  (i)  to  avoid  step  functions  which 
might  disrupt  the  computation,  (ii)  to  model  more  closely  the  gradual 
cooling  of  the  plasma,  by  heat  conductivity,  as  it  approaches  the  foil. 

5-2.  Simulations  with  YRATE. 

YRATE  is  a  calling  program  written  to  provide  CYLINE  with  the 

necessary  input  when  one  wants  to  simulate  plasmas  in  steady-state  or 

with  simple  expansion  characteristics.  It  was  convenient  as  a  debugging 

tool,  prior  to  lengthy  runs  with  LILAC,  to  make  sure  that  (i)  CYLINE 

gave  reasonably  physical  results  for  any  value  of  (T  ,  N  )  within  its 

e  e 

range,  (ii)  some  simplifications  made  to  cut  the  computation  time  had  a 
negligible  effect  on  the  results  (taking  into  account  the  fact  that  this 
is  a  rather  approximate  model  anyway).  For  more  convenience,  the 
program  includes  several  graphic  options. 


5-2-a.  Steady-state  plasma. 
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Let  us  give  some  examples  of  what  the  program  can  do  and  how  it 
models  a  steady-it-.fce  aluminum  plasma.  We  may  make  the  following 
parameter  comparisons: 

1.  Relative  populations  of  the  various  ionic  stages  for  a  given  set 

of  conditions  (T  ,  N  ) .  Figure  5-2  shows  that,  in  the  range  of 
e  e 

densities  relevant  to  this  study,  the  plasma  is  mainly  fully-ionized  at 

1200  eV,  hydrogenic  at  600  eV,  and  beryllium-like  at  60  eV.  It  confirms 

that  the  boundaries  (60  eV<T  <3  keV)  imposed  by  the  availability  of  rate 

e 

coefficients  were  not  restrictive  for  our  purpose.  At  lower 
temperatures,  the  program  still  computes  the  overall  ionic  distribution 
using  Eq.5-1,  but  not  the  populations  of  the  excited  states  of  1-  and 
2-electron  ions. 

2.  Relative  intensity  of  one  line  versus  T  for  a  given  density  N  . 

e  e 

Figure  5-3  shows  the  intensity  of  the  resonance  lines  (a)  La(H-like) , 
(b)  w(He-like),  for  an  optically  thin  plasma.  The  conclusion  is  similar 
to  that  drawn  from  Fig. 5-2,  i.e.  there  is  an  optimal  temperature, 
somewhat  independent  of  density,  for  the  line  emission  from  different 
ionic  stages. 

3.  Ratio  of  corresponding  transitions  in  1-  and  2-electron  ions 

versus  T  for  a  given  density  N  .  Figure  5-4  shows  the  intensity 
e  e 

ratio  of  the  lines  La(H-like)/w(He-like)  between  60  eV  and  1  keV  for 
various  densities  in  the  case  of  (a)  an  optically  thin  plasma,  (b)  an 
optically  thick  plasma  of  depth  <1>»170  ym.  Clearly,  that  ratio 
provides  a  very  sensitive  measurement  of  Tg  (see  Chap. 2)  unless  the 
plasma  density,  hence  the  opacity,  is  high. 

4.  Ratio  of  different  lines  of  the  same  ion  versus  Tg  for  a  given 
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Figure  5-2  :  Ionic  Distribution. 
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Figure  5-3  :  Relative  Line  Intensity  versus  Temperature 


a)  hydrogen-like  resonance  line, 

b)  helium-like  resonance  line. 
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Figure  5-5  :  Intensity  Ratio  ls^-ls  4p  /  ls^-ls  3p. 

a)  optically  thin  plasma, 

b)  optically  thick  plasma. 

The  "inversion  level"  refers  to  the  line  ratio  required 
to  infer  population  inversion  in  the  optically  thin  case. 


^0-8136  305 
UNCLASSIFIED 


DEVELOPMENT  OF  X-RAV  LASER  MEDIA  MEASUREMENT  OF  GAIN 
AND  DEVELOPMENT  OF  C.  .  (U)  ROCHESTER  UNIV  N  V  LAB  FOR 
LASER  ENERGETICS  J  M  FORSVTH  FEB  83 
AFOSR-TR- 83-1136- VOL-1  AFOSR-81-0059  F/G  20/5 


2/2 


F/G  20/5 


density  N@.  Figure  5-5  shows  the  intensity  ratio  of  the  Helium-like 
2  2 

lines  Is  -ls4p/ls  -ls3p  for  (a),  an  optically  thin  plasma,  (b)  an 
optically  thick  plasma.  The  "inversion  level”  refers  to  the  line  ratio 
required  to  infer  population  inversion  in  the  optically  thin  case.  Mote 
that  the  apparent  population  ratio  is  significantly  increased  by 
opacity,  an  important  point  to  keep  in  mind  when  interpreting  the  data. 
The  limit,  in  a  steady-state  plasma,  is  the  blackbody  spectrum  where  the 
emission  from  the  two  lines  would  be  of  approximately  equal  intensity. 


5-2-b.  Expanding  plasma. 

YRATE  can  also  simulate  plasmas  with  simple  expansion 

characteristics,  of  the  form  N  ^Tx,  and  a  cooling  rate  R  *-dT  /dt. 

e  e  c  e 

Depending  on  the  value  of  R  ,  the  atomic  rate  equations  may  not  reach 

c 

their  steady-state  solution  before  the  plasma  conditions  change.  This 

is  exactly  the  situation  we  are  interested  in,  because  a  population 

inversion  may  occur  only  if  the  plasma  cooling  time  is  short  compared  to 

its  recombination  time  (see  Chap. 2).  The  electron  density  is  such  that 

N  /(N  )  =  (z  /z)m  where  z  is  the  radial  coordinate  (see  Figure  5-1), 
e  e  o  o 

(N  )  and  z  the  initial  conditions,  and  m*l,2,3  for  planar,  cylindrical 
e  o  o 

or  spherical  expansion.  If  the  equation  of  state  is  such  that  N'vt*,  it 

e  e 

follows  that  T  /(T  )  »(z  /z)m/x  and  the  cooling  rate  is 
e  e  o  o 


F  *  -(dT  /dz)(dz/dt)  =  (m/x)  (T  /z)  (z  /z)1+^X^v 
c  e  e  o  o  exp 


(5-17-a) 


where  v  is  the  expansion  velocity.  If  (R  )  is  the  initial  cooling 
exp  c  o 


rate, 


R  -  (R  )  [N  /  (N  )  J 


1/m  +  1/x 


(5-17-b) 


As  will  be  seen  in  the  next  section,  typical  numbers  in  the  1-D 

14  2 

simulation  of  point-focus  experiments  on  aluminum  targets  (1=5  10  w/cm 

O 

at  A-0.351  ym)  are  m-3,  z  -200  ym,  (T  )  =2  keV,  v  =10  cm/sec.  For  an 

o  e  o  exp 

adiabatic  expansion  (x«3/2) ,  the  cooling  rate  would  be  R  -20  [N  /(N  )  ] 

c  e  e  o 

eV/psec.  However,  heat  conductivity  maintains  the  temperature  high  in 

the  coronal  plasma,  and  expansion  cooling  is  not  effective  on  very  short 

time  scales  (Forsyth,  et  al,  1976).  Larger  values  of  x  are  required  to 

adequately  describe  the  expansion.  For  x*3,  R^-io  [Ne/(Ne)Q]2//3  eV/psec. 

About  200  ym  above  the  target  surface,  the  density  has  dropped  by  an 

order  of  magnitude,  and  the  cooling  rate  is  of  order  1  eV/psec.  The 

temperature  is  half  its  initial  value.  Taking  the  plasma  at  that  point 
21  20  19  -3 

(T  -1  keV,  N  -10  ,  10  or  10  cm  ) ,  Figure  5-6  shows  the  intensity 

e  e 

2  2 

ratio  of  the  helium-like  lines  Is  -ls4p/ls  -ls3p  for  an  optically  thin 
aluminum  plasma,  with  different  cooling  rates.  If  the  cooling  is  only 
that  due  to  the  spherical  expansion  (5-6-a) ,  population  inversion  occurs 
at  relatively  low  temperature,  i.e.  after  a  large  density  drop, 
especially  if  the  initial  density  is  high.  When  additional  cooling  is 
introduced  (5-6-b) ,  inversion  occurs  at  an  earlier  stage.  The  limit  is 
reached  when  the  cooling  is  so  fast  (5-6-c)  that  the  ionic  populations 
are  virtually  frozen  in  their  initial  distribution. 

5-3.  Simulations  with  LILAC. 

After  having  tested  the  subroutine  CYLINE  (which  contains  our 


computer  modeling)  with  the  simple  YRATE  code,  we  ran  full-scale, 
time-dependent  simulations  with  the  1-dimensional  laser  fusion  code 
LILAC  (op.  cit.).  We  describe  here  the  interface  LILAC-CYLINE,  and 
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in  spherical  expansion*  with  various  cooling  rates  (see  text)* 
to  be  read  from  right  (initial  density)  to  left  (final  density) . 


simulations  of  (i)  point-focus  experiments  on  aluminum  targets  with 
either  IR  or  UV  illumination,  (ii)  line-focus  experiments  on  fluorine 
(UV  illumination  only) .  Results  are  compared  with  the  experimental 
data. 

5-3-a.  Code  LILAC  and  interface  with  CYLINE . 

LILAC  (op.  cit.)  is  a  numerical  code  developed  mainly  to  simulate 
spherical  implosions  of  DT-filled  micro-balloons  for  taertial 
confinement  fusion  (ICF) .  It  computes  the  hydrodynamic  motion  of  not 
only  the  imploding  fuel  and  shell,  but  the  expanding  plasma  as  well.  By 
increasing  the  wall-thickness  of  the  balloons,  it  can  simulate 


laser-matter  interaction  experiments  on  solid  spheres.  In  the  case  of 
point-focus,  the  target  was  represented  by  a  sphere  of  radius  200  ym, 
twice  as  large  as  the  laser  focal  spot  (see  Fig. 5-1),  being  subjected  to 
the  same  ir radiance  as  in  the  experiment  but  with  higher  total  energy. 
The  calculated  X-ray  intensity  was  then  scaled  to  the  measured  laser 
pulse  energy.  This  procedure  is  known  to  give  best  results  when  trying 
to  simulate  with  a  spherically  symmetric  code  the  inherently 
two-dimensional  behavior  of  flat  target  irradiation  (Delettrez  1978-82) . 
The  same  method  was  used  in  the  line-focus  case,  with  a  simple  intensity 
scaling. 

The  code  uses  Lagrangian  fluid  dynamics.  "The  local  density  of  the 
species  is  found  from  the  instantaneous  volume  of  each  cell  in  [a] 
spatial  grid.  (...)  The  logic  of  the  code  (...)  requires  that  all 
species  have  the  same  macroscopic  velocity  and  that  all  of  the  ions  have 
the  same  temperature.  Only  one  momentum  equation  is  then  used,  but  two 
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energy  equations,  one  for  the  electrons  and  one  for  the  ions,  are  still 
necessary.  Energy  exchange  mechanisms  include  flow  work,  thermal 
conductivity,  coulomb  collisions  between  the  electrons  and  ions,  (...) 
and  laser  absorption.  The  thermal  conductivities  are  flux  limited  with 
the  value  of  the  limiter  to  be  specified  by  the  user"  (Goldman  1976) . 
Nultigroup  suprathermal  electron  transport  (Delettrez  &  Goldman  1976) 
and  multigroup  radiation  transport  with  opacities  obtained  from  Local 
Thermodynamic  Equilibrium  tables  (Huebner,  et  al,  1977)  are  also 
included. 

The  code  can  simulate  experiments  with  (i)  IR  light  (1.054  ym) 

where  coupling  of  the  radiation  field  to  the  target  is  due  mainly  to 

resonant  absorption,  with  production  of  suprathermal  electrons,  (ii)  UV 

light  (0.351  ym) ,  where  inverse  bremsstrahlung  is  the  dominant 

absorption  mechanism  (see  Sec. 4-1) .  The  temperature  drop.  At  , 

e 

introduced  in  CYLINE  to  simulate  the  heat  sink  effect  of  the  foil  in 

front  of  the  target  (see  Sec.5-l-c)  translated  into  a  loss  term  in 

LILAC'S  energy  equation  for  the  electrons:  AE»n  C  At  ,  where  n  is  the 

eve  e 

number  of  electrons  in  a  zone,  and  the  heat  capacity  for  electrons. 
This  is  equivalent  to  a  radiation  loss  by  the  plasma.  It  was  on  that 
basis  that  LILAC  computed  a  new  electron  temperature  for  the  next  time 
step. 

The  simulations  were  run  with  a  gaussian  laser  pulse,  fwhm«0.5 
nsec,  peaked  at  t*0.644  nsec  (initial  time:  t«0) .  Depending  on  the 
laser  intensity,  hence  the  plasma  expansion  velocity,  the  simulations 
were  stopped  after  2  to  4  nsec,  enough  time  for  most  of  the  plasma  to  be 
at  distances  21  an  from  the  target  surface,  where  the  density  is  too  low 
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for  relevant  experimental  observations  (at  least  in  our  set-up).  The 

output  listings,  at  time  intervals  0.2  nsec,  included  the  plasma 

parameters  for  each  cell  of  the  spatial  grid  (in  particular  (T  ,N  ), 

©  © 

(T^KL),  v#xp) ,  and  several  features  specific  to  the  subroutine  CYLINE: 

1.  For  each  cell,  the  population  densities  of  bare  nuclei,  and  1-  and 
2-electron  ions  in  the  ground  state  and  excited  states  up  to  n*5  (n, 
principal  quantum  number).  The  average  degree  of  ionization,  <Z>,  was 
also  computed. 

2.  The  cumulative  photon  count  on  a  "detector"  with  10  ym  spatial 
resolution,  placed  at  90°  from  the  laser  axis  (see  Sec.5-l-b) .  The 
count  was  given  for  each  line  of  the  resonance  series  n-+l  (n*2  to  5)  of 
hydrogenic  and  helium-like  ions.  Results  were  given  for  (i)  an 
optically  thin  plasma,  (ii)  an  optically  thick,  Doppler-broadened  plasma 
(see  Sec.5-l-b-2) . 

3.  The  optical  gain,  both  instantaneous  and  time-averaged  (from  the 
beginning  of  the  pulse)  on  the  resonance  transitions  just  mentioned, 
with  similar  spatial  resolution.  The  gain  was  computed  in  the  direction 
perpendicular  to  that  of  the  laser  axis  (along  the  line-focus  if 
applicable) . 

A  graphic  routine  was  written  to  display  and  compare  the  results. 

5-3-b.  Simulation  of  point-focus  experiments  on  aluminum. 

* 

Figure  5-7  shows  the  simulation  of  shot  #2789  (see  Fig. 4-10),  when 
UV  light  (1*5  lO^Sf/cm^)  was  focused  at  the  surface  of  an  aluminum  slab. 
A  100  ya  foil,  placed  200  ym  in  front  of  the  slab,  acted  as  a  heat  sink. 
The  aain  experiaental  features  (increased  X-ray  emission  near  the  heat 
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sink,  anomalous  line  ratios,  see  also  Figure  5-8),  are  well  reproduced 

when  the  "characteristic  cooling  time"  t  is  set  around  30  psec.  It 

c 

corresponds  to  an  enhanced  cooling  rate  R*,  in  addition  to  that  due  to 

the  expanding  plasma  (see  Sec.5-2-b).  In  the  vicinity  of  the  foil, 

R**-dT  /dt^T  /T  .  On  the  side  of  the  foil  nearest  to  the  target,  T  is 
c  e  e  c  e 

of  order  1  keV  and  R*-30  eV/psec.  According  to  Pig. 5-6,  it  is  about  the 

c 

20  -3 

rate  required  to  "freeze"  a  plasma  of  density  N  £  10  cm  and  to  produce 

e 

a  population  inversion  in  the  helium-like  ion. 

As  shown  in  Fig. 5-7,  the  temperature  drop  near  the  heat  sink 

produces  a  pressure  gradient.  The  pressure  exerted  by  the  upstream 

plasma  creates  a  dense  region  downstream  after  the  foil,  appropriate  for 

enhanced  recombination.  The  gain  region  travels  with  the  velocity  of 

8  * 

the  expanding  plasma  (=:10  cm/sec  for  shot  #2789  ) .  While  the  local  gain 

averaged  over  the  plasma  recombination  time  does  not  exceed  10  cm-1,  the 

instantaneous  gain  can  be  much  higher  ('vlOO  cm”1  over  <100  psec).  The 

gain  x  length  product  remains  small  because  of  the  plasma  dimensions,  a 

few  hundred  microns  in  the  transverse  direction.  Fig. 5-8  is  a  side  by 

side  comparison  of  the  simulation  (5-8-a,  same  as  5-7-a)  and  the 

2 

experimental  data  (5-8-b,  similar  to  4-10  limited  to  the  Is  -ls4p  & 

2 

Is  -ls3p  lines,  and  5-8-c,  densitometer  scans) . 

Let  us  add  a  few  comments  about  the  simulation.  Intensity 

2  2 

inversion  between  the  lines  Is  -ls4p  and  Is  -ls3p  is  not  observed 

(i)  when  opacity  is  not  included.  Far  from  the  target,  the  density  is 

low  but  the  size  of  the  plasma  is  large  and  self-absorption  is  still 

2  2 

important.  It  affects  Is  -ls3p  more  than  Is  -ls4p,  m  ratio  of  the 


lines'  oscillator  strengths  (0.153  vs.  0.057).  Opacity  affects  most 


89 


strongly  the  dense  region  near  the  target  surface.  If  it  were  not  taken 

into  account,  the  ratio  of  the  two  peaks  of  the  emission  (i.e.  at  the 

target  surface  and  near  the  heat  sink)  would  be  much  higher  than 

observed  experimentally,  by  almost  an  order  of  magnitude. 

(ii)  when  there  is  no  heat  sink,  i.e.  r  ■«.  The  line  intensities 

c 

decrease  steadily  with  distance.  During  the  expansion,  their  ratio 

remains  approximately  what  it  was  at  distance  100  ym  from  the  target. 

The  consequences  of  a  different  location  for  the  heat  sink  are  apparent 

on  Fig. 5-6.  If  the  foil  is  closer  to  the  target  and  interacts  with  a 

denser  plasma,  a  very  high  cooling  rate  is  required  to  "freeze”  the 

ionic  distribution.  If  the  foil  is  farther  away,  intensity  inversion 

may  still  occur  but  at  low  density,  hence  low  gain.  Furthermore  these 

geometries  are  somewhat  impractical  experimentally  (see  Sec. 4-2). 

There  is  a  large  margin,  factor  ^2-3,  for  the  choice  of  t  . 

c 

However,  if  T  is  too  long,  the  cooling  is  not  efficient  enough  to 
c 

produce  an  inversion.  If  T  is  too  short,  not  only  may  it  become 

c 

unphysical  as  a  heat  transfer  time  constant,  but  on  practical  grounds 

the  computation  cannot  proceed  properly:  the  temperature  Tg  quickly 

drops  to  the  level  where  available  rate  equations  do  not  apply  (see 

Sec.5-l-a)  and  the  photon  count  for  the  transitions  of  interest  has  to 

be  temporarily  bypassed.  In  our  simulations,  we  used  t  as  a  free 

c 

parameter  to  get  an  approximate  match  with  the  experimental  data.  While 

the  cooling  rate  R*  contains  a  space  and  temperature  dependence,  it 
c 

would  certainly  be  an  improvement  if  LILAC  were  to  do  the  calculation 
itself,  taking  into  account  all  plasma  parameters  at  each  time  step.  It 
would  have  to  calculate  the  behaviour  of  a  hot  plasma  expanding  toward  a 


cold  surface.  That  may  be  possible  in  the  near  future,  due  to  the 
current  interest  in  double-shell  targets  for  laser  fusion.  In  our  case, 
the  incident  laser  would  "ignore"  the  outer  shell  and  create  a  plasma  on 
the  surface  of  the  inner  shell.  The  outer  shell  would  then  act  as  a 
heat  sink  for  the  expanding  plasma. 

LILAC  was  also  run  in  its  IR  version  (X=1.054  ym) .  Figure  5-! 

shows  the  simulation  of  shot  #972  (see  Fig. 4-5)  when  IR  light 
15  2 

(1=10  W/cm  )  was  focused  at  the  surface  of  an  aluminum  target  simi  s 

* 

to  that  used  later  in  shot  #2789  .  As  mentioned  in  Chap. 4,  when  IR 

instead  of  UV  light  is  used,  (i)  the  density  of  the  plasma  formed  at 

21 

the  surface  of  the  target  is  one  order  of  magnitude  lower  (10 
22  -3 

vs.  10  cm  ),  (ii)  absorption  of  the  light  being  less  efficient,  higher 

incident  laser  intensities  are  required  to  produce  a  population 

inversion  in  the  helium-like  ion.  At  comparable  levels  of  absorbed 

intensity,  the  smaller  number  of  absorbers  leads  to  a  somewhat  higher 

temperature  in  the  IR  case  (see  Fig.5-9-b  vs.  5-7-b) .  The  combination 

of  lower  density  and  higher  temperature  reduces  the  rate  of  three-body 

recombination  and  makes  it  more  difficult  to  achieve  significant  gain. 

For  this  IR  simulation,  the  "characteristic  cooling  time"  t  was 

c 

set  at  the  same  value,  30  psec,  as  previously  (UV  case) .  Basic 
experimental  features  are  reproduced  (increased  X-ray  emission  near  the 
heat  sink,  enhanced  line  ratios).  While  population  inversion  does  occur 
in  the  expansion,  there  is  no  intensity  inversion,  one  main  reason  being 
the  absence  of  opacity,  due  to  the  lower  density.  Another  consequence 
of  lcwer  density  is  the  smaller  pressure  gradient  produced  by  the 
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Figure  5-9  :  Simulation  of  Shot  #972  ( I R>  Point-Focus) 


temperature  drop  near  the  heat  sink.  The  dense  region  that  it  creates 
downstream  after  the  foil  (see  UV  case,  Fig.5-7-b)  is  barely  visible  in 
the  IR  case  (Fig.5-9-b). 


Finally,  we  should  point  out  that,  in  both  the  UV  and  IR 
simulations,  the  X-ray  yield  and  its  distribution  among  the  spectral 
lines  of  1-  and  2-electron  ions  were  in  general  agreement  with  the 
experimental  data  (factor  ^2  for  the  absolute  yields,  and  $30%  for  the 
line  ratios) .  The  fit  is  reasonable,  considering  the  approximations  in 
the  model  and  the  uncertainty  in  the  film  calibration. 

5-3-c.  Simulation  of  line-focus  experiments  on  teflon. 

Some  changes  were  required  on  LILAC  and  CYLINE  to  run  simulations 
of  shots  on  teflon. 

(i)  Since  the  table  of  materials  included  in  LILAC  does  not  contain 

data  for  fluorine,  we  used  instead  of  (CF„)  a  fictitious  CO_  target  at 

2  n  2 

solid  density  and  room  temperature  (actually  0°C) .  The  similarity  in 
composition  and  atomic  number  (Z*8  for  0,  9  for  F)  should  not  lead  to 
very  different  results. 

(ii)  A  separate  version  of  CYLINE  (file  name  CYOXYG)  was  used,  with 
atomic  data  for  oxygen  and  the  proper  scaling  laws. 

* 

Figure  5-10  shows  the  simulation  of  shot  #4471  (see  Fig. 4-15), 

13  2 

when  UV  light  (1=1.2  10  W/cm  )  was  focused  on  a  100  ym  x  2  mm  area  of 
teflon  (i.e.  "CO^").  A  100  ym  foil,  placed  200  ym  in  front  of  the  slab 
acted  as  a  heat  sink.  Fig.5-10-a  compares  the  emission  on  the 
hydrogenic  transitions  L. (ls-3p)  and  L  (ls-2p) .  The  "characteristic 

p  a 

cooling  time"  was  set  at  300  psec.  There  is  less  cooling  in  the 


vw. 
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line-focus  geometry,  compared  to  point-focus,  because  the  plasma  has  a 
smaller  volume  of  nearby  vacuum  in  which  to  expand.  A  periodic  spatial 
variation  in  the  pumping  radiation  should  help  solve  that  problem 
(Silfvast  et  al  1979-a) . 

Here  again  the  main  experimental  features  are  well  reproduced. 

X-ray  emission  increases  near  the  foil,  and  is  even  stronger  than  at  the 


target  surface  (in  agreement  with  the  pattern  recorded  for  L  ) ,  due  to 

3 

lower  opacity  and  lower  temperature.  There  is  less  opacity  near  the 


foil  because  the  relative  change  in  plasma  length  is  small  compared  to 


the  drop  in  density.  The  lower  temperature  is  more  adequate  for 

bare-nuclei  recombination  than  the  hot  focal  region.  The  intensity 

inversion  (L<LQ)  in  the  expansion  is  the  combined  result  of  (i)  a 
a  p 

population  inversion,  (ii)  a  difference  in  opacity  (the  oscillator 
strengths  for  1^  and  are  in  the  ratio  0.42  vs.  0.08).  The  gain 

region,  while  moving  in  time,  is  more  stable  than  in  the  simulation  of 

* 

shot  #2789  (Fig. 5-7:  point-focus,  aluminum  target)  because  (i)  the 

7  8 

expansion  velocity  is  lower  (3x10  vs.  10  cm/sec),  due  to  the  lower 


laser  intensity,  (ii)  the  population  inversion  occurs  earlier  and  is 


less  dependent  on  the  propagation  of  a  high  density,  low  temperature 
plasma . 


The  line  intensities,  however,  decrease  much  faster  than  it  appears 
in  the  data.  This  points  to  the  difficulty  of  describing  with  a 
1-dimensional  code  a  line-focus  experiment  where  the  plasma  expansion  is 
cylindrical  (or  possibly  planar,  see  Fig.4-15-a),  rather  than  spherical. 
Our  pseudo-planar  geometry  using  a  large  sphere  (Fig. 5-1)  can  model  the 
plasma  overall  divergence  but  not  its  asymetry  and  the  associated  (known 


and  unknown)  hydrodynamic  effects.  In  the  vicinity  of  the  foil,  the 
gain  averaged  over  the  plasma  recombination  time  is  about  0.1  cm“*.  its 
instantaneous  value  can  reach  5  cm*' 1  over  ^  200  psec  for  a  gain  x  length 
products  1. 

In  this  chapter,  we  have  tried  to  model  the  experimental 
observation  of  intensity  inversions  in  laser-produced  plasmas  (IR  and  UV 
illumination,  point-focus  and  line-focus) .  Using  existing  laser  fusion 
code  (Delettrez  1978-82)  and  rate  equations  (MCWhirter  6  Hearn  1963) ,  we 
found  good  agreement  with  the  data  when  we  included  (i)  enhanced 
recombination,  via  cooling,  (ii)  plasma  opacity.  The  gain  region 
travels  with  the  velocity  of  the  expanding  plasma.  Its  value,  high 
locally  in  time  and  space,  seems  below  the  requirements  for  a 
time-integrated  experimental  observation. 
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Summary 

This  thesis  describes  the  progress  made  in  the  development  of  an 
XUV  amplifier  based  on  the  recombination  of  a  laser-produced  plasma. 
Experimental  results  and  numerical  simulations  have  led  to  a  better 
understanding  of  the  process  leading  to  a  population  inversion  during 
the  plasma  expansion.  While  several  important  questions  require  further 
study,  direct  demonstration  of  coherent  amplification  at  'V100  &  seems 
possible  in  the  forseeable  future.  We  review  our  results,  their 
limitations,  and  the  prospect  for  further  study. 

Population  inversions  were  observed  in  the  helium-like  ions  Al  & 

+10  |  Q 

Mg  and  the  hydrogenic  ion  P  ,  by  intoducing  a  foil  in  the  plasma 
expansion.  While  the  idea  in  itself  is  not  new  (Bhagavatula  &  Yaakobi 
1978),  (i)  we  demonstrated  that  the  inversion  was  caused  neither  by 

resonance  mechanisms  such  as  photopumping  or  charge  transfer,  nor  by 
channelling  of  the  plasma  flow,  (ii)  we  defined,  after  a  parametric 
study,  the  optimum  target  design  (a  slit-foil  of  an  arbitrary  material, 
200  pm  wide  and  200  ym  from  the  original  plasma  formation) ,  and 
illumination  conditions  (short  wavelength  laser,  X»0.351  ym,  of 
intensity  in  the  range  1013-1015  Watts/cm2  for  targets  of  atomic  number 
9  to  13).  Our  experiments  used  a  pump  laser  focussed  on  the  target  in 
either  a  point-focus  or  a  line-focus  (using  specially  manufactured 
cylindrical  optics).  Inversion  was  observed  in  both  cases,  with 

|  appropriate  target  materials  and  laser  intensities. 

I 

i 
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Simulations  of  the  experiments  were  done  by  combining  existing 
1-dimens ional  hydrodynamic  code  (Delettrez  1978-82)  6  rate  equations 
(McMhirter  k  Hearn  1963)  with  our  model,  which  included  (i)  enhanced 
recombination,  via  cooling  of  the  plasma  in  the  vicinity  of  the  foil, 
(ii)  radiation  transport  (opacity).  The  main  features  of  the 
time-integrated  data  (overall  X-ray  flux,  increased  X-ray  emission  near 
the  foil,  anomalous  line  ratios)  could  be  reproduced  with  the  proper 
choice  of  a  "cooling  rate”.  Lowering  the  temperature  of  a  dense, 
highly-ionized  plasma  results  in  a  recombination  cascade  dominated  by 
collisions,  which  tend  to  populate  high-lying  quantum  states  of  the 
recombining  ion  (Gudzenko  &  Shelepin  1965) .  The  simulations  predict 

that  the  gain  region  travels  with  the  velocity  of  the  expanding  plasma. 

2  -1 

The  gain  could  be  high  locally  in  time  and  space,  10  cm  over  £  100 
psec  on  Al  (3  D-4  F)  at  129.7  A  and  5  cm  over  ^200  psec  on 
F  (Balmer  a)  at  80.9  with  point-focus  and  line-focus  illumination 
respectively.  However,  the  transverse  dimensions  of  the  plasma  (about 
200  ym  and  2.5  mm  in  the  region  of  highest  gain)  were  such  that  the  gain 
x  length  product  was  approximately  1  in  both  cases.  This  is  3-5  times 
less  than  what  is  required  for  a  direct  gain  measurement  from  recordings 
of  the  amplified  spontaneous  emission  (Slaymaker  1978). 

Our  study  was  limited  in  several  respects:  (i)  the  data  were 
time-integrated ,  (li)  the  simulations  were  done  with  a  1-dimensional 
code,  (iii)  a  phenomenolog i c  description  was  used  for  the  interaction  of 
the  expanding  plasma  with  the  foil.  This  suggests  directions  for 


further  research 
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Measurements  with  an  X-ray  streak-camera  (Letzering  1980)  would 
give  time-resolved  information  about  the  X-ray  emission  near  the  foil 
and  the  population  inversion  in  the  expansion.  It  would  become  possible 
to  verify  some  of  the  code  predictions  such  as  that  of  a  traveling 
inversion  zone.  Direct  gain  measurements  are  not  in  order  at  this  point 
on  a  time-integrated  basis  (gain  x  length  products  are  <0.3,  averaged 
over  the  plasma  lifetime).  With  the  time  resolution,  they  could  become 
feasible,  considering  that  (i)  we  used  only  a  fraction  of  the  pump 
intensity  now  available  at  the  laboratory  (20  J  vs.  60  J  of  uv  light  in 
500  psec) ,  (ii)  further  upgrade  of  the  pump  laser  is  forthcoming,  in 
both  single-beam  (GDL)  and  multi-beam  (frequency-tripled  OMEGA) 
configurations,  to  almost  250  J  in  1  nsec,  (iii)  progress  in  target 
design,  hence  gain,  should  be  expected.  With  larger  lasers,  one  could 

not  only  use  a  longer  line-focus  to  increase  the  gain  x  length  product, 
but  also  investigate  higher-Z  targets  in  order  to  scale  the 
recombination  scheme  toward  shorter  (i.e.  X-ray)  wavelengths. 

Preliminary  experiments  with  titanium  (Z*22)  were  encouraging,  showing 

2  2  +20 
an  anomalous  line  ratio  Is  -ls3p/ls  -ls4p  in  the  recombining  ion  Ti 

The  modeling  could  be  improved  in  a  variety  of  ways: 

(i)  More  accurate  rate  coefficients.  Atomic  constants  are  difficult 
to  calculate  except  for  the  hydrogenic  ion,  but  updated  values  supported 
by  experimental  data  are  regularly  proposed  (see  for  example  Landshoff  & 
Perez  1976) . 

(ii)  Multi-dimensional  hydrodynamic  code.  Since  our  target  design 
does  not  have  spherical  symmetry,  1-dimensional  modeling  is  only 
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approximate.  New  codes  developed  for  inertial  confinement  fusion 
research,  such  as  the  2-dimensional  codes  LASNEX  (Zimmerman  1973) ,  used 
by  Hagelstein  (1981),  or  SAGE  (Craxton  1978-82),  would  permit  a  more 
detailed  description  of  the  time-  and  space-dependent  plasma  parameters. 

(iii)  Plasma  cooling.  Our  phenomenologic  approach  could  be  replaced 
by  a  more  rigorous  treatment  combining  plasma  physics,  thermodynamics, 
and  numerical  analysis.  It  was  beyond  the  scope  of  this  thesis  to 
calculate  the  behaviour  of  a  hot  plasma  expanding  toward  a  cold  surface, 
a  problem  as  complex  as  the  absorption  of  light  in  laser-matter 
interaction  experiments,  still  only  partially  understood  despite 
extensive  studies  in  the  laser  fusion  community  (Godwin  1979,  Max  & 
Estabrook  1979) . 

In  summary,  we  present  in  this  thesis  strong  evidence  of  population 
inversion  in  a  laser-produced  plasma,  due  to  enhanced  recombination  via 
cooling.  That  conclusion  is  drawn  not  from  an  isolated  observation,  but 
from  an  extended  set  of  experimental  data,  whose  main  features  are  in 
general  agreement  with  the  predictions  of  an  advanced  hydrodynamic  code. 

Gain  measurements  seem  possible  in  the  near  future  because  the  two  main 
problems,  size  of  the  pump  laser  and  time  resolution,  can  be  addressed 
with  present  or  forthcoming  technologies.  The  projected  energy 
available  at  0.351  pm  from  either  GDL  or  OMEGA  (6  beams)  being  a  factor 
.10  higher  than  what  we  used  in  this  work,  a  1-2  cm  long  line-focus 
should  be  possible  instead  of  2  mm.  The  travel  time  of  a  photon  in  such 
a  medium,  30  psec/cm,  is  much  shorter  than  the  duration  of  the  inversion 
(200  psec  in  F  ).  Current  X-ray  streak  cameras  (op.cit.)  are  already 
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faster  than  required  by  two  orders  of  magnitude  (^1  psec) .  The 
increased  X-ray  intensity  in  a  scaled-up  experiment  should  compensate 
for  the  lower  level  of  exposure  in  a  time-resolving  detection  system. 

Let  us  sketch  an  experiment  which  could  demonstrate  gain  on  F4® 

(Balmer  a )  at  80.9  &.  A  serrated  target  would  be  manufactured,  made  of 

stripes  of  teflon  separated  by  deep  grooves.  Stripes  and  grooves  would 

each  be  M  mm  wide.  200  Joules  of  UV  light  (0.351  pm),  delivered  in  1 

nsec,  would  be  focused  on  that  target  over  a  line  (2  cm  x  100  Pm) 

transverse  to  the  direction  of  the  stripes.  The  incident  laser,  of 
13  2 

intensity  10  W/cm  ,  would  produce  ten  individual  fluorine  plasmas  (1 

mm  x  100  pm) .  Two  100  pm  sheets  of  titanium  would  be  set  to  form  a  slit 
200  pm  wide  and  200  un  above  the  teflon,  in  a  direction  parallel  to  that 
of  the  line-focus.  Such  a  design  would  allow  lateral  cooling  of  the 
plasma,  through  expansion  or  contact  with  a  heat  sink,  while  maintaining 
the  line-geometry  required  for  an  ASE  diagnostic.  Two  grating 
spectrographs  would  record  the  80.9  A  radiation  emitted  along  the 
line-focus  and  in  the  transverse  direction.  The  film  plane  of  each 
spectrograph  (Fig. 3-5)  would  be  replaced  by  the  entrance  slit  of  an  xuv 
streak  camera.  The  experiment  would  be  repeated  with  a  90°  rotation  of 
both  the  line-focus  and  the  target  in  order  to  cross-reference  the  two 
spectrographs.  Additional  data  could  be  taken  at  other  azimuths.  The 
gain  x  length  product  is  expected  to  be  large  enough  (  ^  5  over  200  psec) 
to  be  detected  by  means  of  the  angular  variations  of  the  XUV  emission. 
Such  project  would  require  a  significant  engineering  effort,  but  this 
thesis  has  shown  that  the  fundamental  obstacles  have  been  overcome  and 


that  a  direct  demonstration  of  XUV  amplif ication  could  be  achieved  soon. 
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Appendix  A 


Photometry  and  film  calibration 


Photometry  can  be  an  accurate  science,  but  is  especially  difficult 
when  it  involves  use  of  film.  Considerable  effort  was  made  to  define  a 
standard  procedure  to  reduce  our  data,  making  it  reasonable  to  compare 
the  results  of  different  experiments.  This  is  important  for  the  crystal 
spectrograph  because  the  set-up,  adapted  to  the  particular  goal  of  each 
experiment,  can  affect  the  film  exposure  in  a  number  of  ways:  dimension 
and  position  of  the  source,  width  of  the  resolving  slit,  angle  of 
observation,  types  of  filter,  crystal,  film,  etc...  It  is  also 
necessary  to  be  consistent  in  the  way  the  microdensitometry  is  done. 

This  appendix  presents  the  procedure  we  used,  our  film  calibration  and  a 
comparison  with  pubished  data.  A  computer  program  (filename  PHOTON)  was 
written  to  handle  that  information.  At  this  stage  of  the  X-ray 
amplifier  project,  the  grating  spectrograph  was  used  more  to  obtain 


qualitative  than  quantitative  results.  We  present  its  imaging 
characteristics  but  no  film  calibration. 


A-l.  Data  reduction  procedure. 

First  one  must  establish  the  film  exposure  geometry.  With  the 

coordinates  (x,y,z)  and  (y',z')  indicated  on  Figure  3-4  and  M  *z'/z  the 

z 

transverse  magnification  between  the  source  and  the  film,  the  time 
dependent  population  density  n^Ujt)  of  the  upper  level  of  the  recorded 
transition  is  related  to  the  time  integrated  illumination  E 
(photons/unit  area)  on  the  film  as  follows: 


"  I.v'.v 
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line 


5tis  the  integrated  reflectivity  of  the  crystal ,  £  the  transmittance  of 
the  filter (s) ,  w  the  width  of  the  slit,  L  the  distance  traveled  by  the 
X-rays  from  the  source  to  the  film,  A  the  Einstein  coefficient  for  the 
transition,  and  k  the  opacity  factor  in  the  direction  of  observation. 
Except  for  w,  these  quantities  are  wavelength  dependent.  Eq. (A-l) 
assumes  that  the  acceptance  angle  in  the  spectral  direction  is  crystal 
limited.  In  the  two  extreme  cases,  no  slit  (w  »  apparent  source  size 
£*)  and  narrow  slit  (w  <<  ) ,  the  spatial  integration  over  the  plasma 

reduces  to 


/  dr  n  (r,t)  k(£,t)  w  >>  $>*  (A-l-a) 

u  s 

w  (1+1/M  )  /  dx  dy  n  (r,t)  k(r,t)  .  5  (z-z'/M  )  w  «  0*  (A-l-b) 

z  u  —  —  z  s 

When  a  shadowgram  is  cast  by  the  outer  edge  of  the  slit  or  by  a 

semi-infinite  aperture,  one  gets 

/  dz  /  dx  dy  n  (r,t)  k(£,t)  (A-l-c) 

Z'/Mz 


The  right-hand  side  of  Eq. (A-l)  is  equivalent  to  the  formula  used 
by  Henke  et  al  (1978)  to  compute  the  total  number  N^,  of  photons  emitted 
by  the  source  into  4tt.  That  formula  is  valid  only  if  (i)  there  is  no 
slit,  (ii)  the  emission  from  the  source  is  isotropic,  that  is  to  say 


opacity  does  not  depend  on  the  direction  of  observation.  In  the  most 


general  case  one  can  compute  the  apparent  value  of  in  the 
direction  of  observation,  by  adding  up  the  photons  recorded  in  each 
resolution  element  on  the  film: 

N*  =  (47TL/31.t)  I  f  dy'  E[y'  ,  z'=  jw(M  +1)  ]  (A-2) 

j=0,°°  spectral  z 

line 

Figure  A-l-a  shows  how  the  microdensitometry  is  done.  A 

Joyce-Loebl  model  3C-S  has  been  used.  1  and  1  are  the  dimensions  of 

t  s 

the  area  scanning  the  film  in  the  spectral  direction: 

2 
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I*,  and  L  are  the  dimensions  of  the  adjustable  slit  in  front  of  the 

w  S 

Joyce-Loebl 's  detection  system,  Mq  is  the  magnification  of  the 
microscope  objective  collecting  the  light  (0) ,  "2.2"  comes  from  the 
built-in  choice  of  optical  conjugation  between  film  and  detector.  A 
typical  line  has  the  profile  indicated  on  the  graph  D  versus  Y'.  we 
assume  that  it  can  be  described  as  follows: 

D(Y',Z')  «  Dq  +  (Dm-DQ)  exp[-(2  Y*/ FWHM)2  ln2]  (A-4) 

Since  the  relationship  D(yr  ,z')  ■*-*■  e  (y'  ,z')  between  density  and  exposure 
(number  of  photons/unit  area  of  film)  is  non-linear,  the  most  convenient 
way  to  proceed  is  to  sample  Eq. (A-4)  (see  Figure  A-l-b) .  As  interval  of 
sampling  on  the  film  one  takes  the  width  of  the  scanning  area  (i.e.  ml 

s 

on  the  paper  trace)  and  one  stops  at  y'<fwhm  where  the  value  of  the 
exponential  has  dropped  to  .06: 

pQ  ■  Integer  (fwhm/1  )  *  Integer  [2. 2  (FWHM/L  )  (M  /  m) )  (A-5) 


Density 


Paper 

trace 


For  a  reasonable  sampling  accuracy,  Lg  must  be  chosen  such  that  pQ>5. 

is  chosen  large  enough  to  reduce  the  noise  from  film  granularity 
without  affecting  the  resolution.  It  follows: 

f  dy'  E(y',z')  *  1  Z  {E(Y’=pml  ,Z')  -  E(D=D  )}  (A_6) 

spectral  s  p=-p  p  s  0 

line  0  0 

Let  us  combine  Eq. (A-2)  and  (A-6): 

N*  a  (4TTL/9lf)  1  Z  Z  {  E (Y’=p  m 1  ,  Z  *  ■  j  m  w  (M  +1) )  -  E(D=D  )}  (A-7) 

J-°«‘  3  2  0 

With  an  adequate  film  calibration,  a  digitized  reading  of  optical 
densities  on  the  paper  trace  leads  directly  to  N*. 

A-2.  Film  calibration. 

In  the  crystal  spectrograph  we  used  the  films  Kodak  RAR-2491  (for 
A>4  £)  and  Kodak  No-screen  (A<8  J() .  Before  drawing  the  correspondence 
between  exposure  and  density  for  those  films  at  various  wavelengths,  it 
is  necessary  to  define  what  is  actually  measured  with  the 
microdensitometer.  If  all  the  light  emerging  from  the  film,  i.e. 
direct  and  scattered  components,  are  collected,  one  measures  the 
so-called  diffuse  density.  As  the  numerical  aperture  of  the  objective 
(0)  decreases  (see  Fig.A-l-a),  less  scattered  light  is  collected  and  the 
measured  density  increases.  The  correction,  which  depends  on  the 
particular  scattering  properties  of  each  film,  has  been  modeled  and 
measured  for  No-screen  by  Brown  (1976),  Dozier  (1976)  et  al.  We  checked 
that  correction  and  extended  it  to  RAR-2491.  The  Joyce-Loebl 
densitometer  was  equipped  with  four  different  objectives  (numerical 
apertures  .10,  .25,  .40,  .85)  to  scan  exactly  the  same  region  of  samples 
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from  each  film.  The  results  are  shown  on  Figure  A-2.  In  both  cases  the 

baseline  is  the  same  for  the  four  traces.  Note  for  example  that  on 

+11  2 

No-screen  the  density  at  the  maximum  of  the  Al  (Is  -ls2p)  resonance 
line  ("w")  drops  from  2.1  to  1.4S  when  the  numerical  aperture  of  the 
objective  collecting  the  light  is  changed  from  .10  to  .85.  Our  data  for 
No-screen  fit  reasonably  well  the  correction  "measured  density/diffuse 
density"  suggested  by  Dozier  et  al.  The  correction  required  for 
RAR-2491  is  very  similar  to  the  one  used  for  No-screen. 

The  film  calibration  itself  was  done  by  recording  spectra  through 
Aluminum  filters  of  various  thicknesses  (see  Table  A-l) .  Data  for 
No-screen  were  taken  around  three  wavelengths  and  compared  to  those 
given  by  Dozier  et  al.  Data  for  RAR-2491  were  taken  around  7  %-Q  ^  and 
compared  with  the  calibration  given  for  RAR-2490  by  Benjamin  et  al 
(1977).  The  results  for  No-screen  are  presented  on  Figure  A-3-a.  Our 
H-D  curves,  arbitrarily  positioned  along  the  exposure  scale  to  be  in 
coincidence  with  Dozier's  curves  at  low  optical  density,  have 
significantly  lower  slopes.  Much  closer  agreement  is  found  between  our 
data  on  RAR-2491  and  Benjamin's  data. 

Note  that  the  Aluminum  spectrum  was  recorded  simultaneously  on 
No-screen  and  RAR-2491  (shots  12138  and  2201) .  Fig. A-2  shows  that  the 
line  Al*^(w)  has  about  the  same  density  on  both  films,  suggesting  that 


the  H-D  curves  should  be  close  to  each  other  around  D..,r  =1.5.  This 

diffuse 

is  confirmed  by  other  data  such  as  shots  #2903  and  2904  where  the 


spatially  resolved  spectra  from  identical  step  targets  were  successively 


recorded  on  RAR-2491  and  No-screen.  Possible  causes  of  error,  such  as 


differences  in  filter  thickness  and  crystal  reflectivity  between  the 


Calibration  Experiments. 


Comparison  No-Screen 
vs.  RAR-2491  at  7.15  X: 


Figure  A-3  :  Film  Calibration  with  GDL. 

a-  Calibration  of  No-screen  Film, 
b-  Comparison  of  No-screen  and  RAR-2491  Films 
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spectrographs  used  to  record  the  two  spectra,  were  checked  through 
control  experiments.  One  has  to  take  into  account  the  fact  that 
RAR-2491  is  somewhat  more  sensitive  than  RAR-2490  in  the  blue.  It  may 
very  well  hold  in  the  X-ray  regime  too.  A  translation  of  the  RAR-2490 
H-D  curve  (factor  v  1.5  on  the  logarithmic  scale  for  the  exposures), 
combined  with  the  correction  previously  mentioned  for  No-screen,  make 
the  data  compatible.  That  double  correction  is  shown  on  Figure  A-3-b. 

A-3.  Code  PHOTON. 

The  code  carries  out  the  procedure  described  in  Sec.A-1.  It 

contains  film  calibration  curves  for  No-screen  (1.6  &  to  11  &)  and 

RAR-2490  (1.6  &  to  44  &)  from  Dozier  and  from  Benjamin  as  well  as  our 

own  curves  over  the  same  spectral  range,  using  corrections  consistent 

with  those  described  in  Sec.A-2.  Tables  (Dozier  et  al  1976)  necessary 

for  the  conversion  "measured  density/diffuse  density"  are  also  included. 

That  information  is  used  as  a  grid:  the  exposure  corresponding  to  any 

density  for  any  numerical  aperture  at  any  wavelength  within  the  grid  is 

determined  by  linear  interpolation  from  the  nearest  available  curves. 

The  code  requires  as  input  the  user'  selection  of  a  film 

calibration  and  the  relevant  parameters  of  the  microdensitometer: 

magnification  (M  }  and  numerical  aperture  of  the  objective  collecting 
o 

the  light,  width  (L  )  of  the  slit  in  front  of  the  detector,  scale  of  the 

s 

density  wedge,  magnification  (m)  paper  trace/film.  It  also  requires  for 
each  spectral  line  the  wavelength,  width  (FWHM) ,  height  at  the  peak  and 
at  the  base  of  the  profile  on  the  paper  trace.  The  output  gives  for 
each  line  the  value  of  (5l^/L).NT,  see  Eq.A-2.  $1  and  ^  are  readily 


available  in  the  literature,  L  is  provided  by  another  code  (CRYSPEC,  see 
Appendix  B) . 

A-4.  Photometry  of  the  grating  spectrograph. 

Eq.A-1  given  for  the  crystal  spectrograph  applies  also  to  the 
grating  spectrograph  (see  Fig. 3-5)  if  the  quantity  L/%  on  the  right-hand 
side  is  replaced  by 

d  [  +  r(sina  +  sing) ]  /£  s 

0  0 

£  is  the  efficiency  of  the  grating  in  the  relevant  order  of  diffraction, 
r  its  radius  of  curvature,  a  and  B  the  incidence  and  diffraction  angles 
measured  from  the  surface  of  the  grating,  s  and  dQ  the  width  of  the 
entrance  slit  and  its  distance  to  the  source,  respectively.  The  same 
substitution  can  be  made  in  the  other  equations  of  Sec.A-1.  However  the 
procedure  for  microdensitometry  loses  its  relevance  because  the  lines  on 
the  film  are  much  narrower  than  they  are  in  the  non-focusing  crystal 
spectrograph. 

The  film  used,  Kodak  101-01,  was  not  calibrated  because  the 
diffraction  efficiencies  of  our  gratings  at  different  wavelengths  and 
for  different  orders  were  not  known.  References  to  101-01  in  the 
literature  (Burton  et  al  1973,  Koppel  1974)  are  limited  to  wavelengths 
outside  the  spectral  range  we  are  interested  in  (50  ft-200  R)  and  it  is 
unclear  how  an  extrapolation  should  be  done. 


iridi 
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Appen< 

Experimental  waveli 

Most  of  the  data  presented  in  this  thesis  were  recorded  with 
spectrographs  whose  dispersive  element  was  either  a  flat  (non-focusing) 
crystal  or  a  concave  (focusing)  grating.  When  a  crystal  is  used,  the 
wavelength  scale  depends  on  the  relative  position  of  the  source  with 
respect  to  the  spectrograph;  the  wavelength  determination  has  to  rely  on 
some  known  lines  in  the  recorded  spectrum.  When  a  grating  is  used,  the 
wavelength  scale  is  independent  of  the  source  position  but  needs  to  be 
determined  with  accuracy  to  get  full  benefit  from  the  resolving  power  of 
the  instrument. 

This  appendix  describes  for  each  case  a  method  to  determine  the 
wavelength  scale  from  the  recorded  spectrum  itself  provided  that  the 
wavelengths  of  a  minimum  number  of  lines  are  known.  For  completeness  a 
method  applying  to  the  curved-crystal  spectrograph  is  also  presented. 

The  curvature  ay  be  accidental  (a  "flat"  crystal  improperly 
characterized  or  subject  to  stress  in  its  mount)  or  deliberate,  to 
increase  the  spectral  range  of  the  spectrograph  when  the  crystal  cannot 
be  positioned  near  the  source. 


B-l.  Flat-crystal  spectrograph. 


B-l-a.  Description  of  the  problem. 


Figure  B-l  :  Geometry  of  a  Flat-Crystal  Spectrograph 
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planes  ($*) .  The  normals  to  3  and®  are  assumed  to  be  coplanar.  The 
X-rays  are  diffracted  according  to  Bragg's  law  X*2D.sin(0+i[i)  *2D.sin0*, 
where  D  is  the  interplanar  spacing  of  the  crystal.  One  gets  from 
Fig.B-l-a 

y  =  a  +  b/ [tana  +  tan  (9*+^)  3  +  h  sin  /  sin  (d*-ip)  sin  (a+9*+i/>)  (B-l) 

Positions  y  on  the  film  are  measured  with  respect  to  an  arbitrary  point 
0  whose  distance  to  the  point  Q  is  a.  The  distance  (STS')  that  the  ray 
of  wavelength  X  travels  from  the  source  to  the  film  is 

L=  [b  cos2a  +  2  h  sini|j  cos  (a+0*)/ sin(0*-i|O  3  /  sin(a+0*+^)  (B-2) 

This  distance  is  useful  for  (i)  photometric  calculations,  (ii)  the 
determination  of  the  wavelength-dependent  linear  magnification  of  a  slit 
in  a  spatially-resolving  arrangement. 

If  a  densitometer  of  magnification  m  is  used,  Eq. (B-l)  still 
applies  provided  one  makes  the  following  substitutions:  Y=my,  A=ma, 

B*mb,  H=mh.  In  the  most  general  case  Eq. (B-l)  involves  five  independent 
parameters  (A,  B,  a,  h,  i(i).  The  sensitivity  of  Y  and  L  to  these 
parameters  is  measured  by  the  partial  dervatives: 

9y/9a  =  1 

9Y/9B  =  (tana  + tan (9*+^) 3  1 

3Y/3a  =  -B  (9Y/9B)2/  cos2a  -  h(9y/9h)  /tan(a+9*+^)  (b-3) 

9y/9h  =  sin  2ip  /  sin  (9*-^)  sin(a+9*+^) 

9Y/9iJ>  =  -BOY/9B)2/ cos2  (9 *+![;) 

+H(9y/3H)  [cosa  -  cos  2ip  cos(a+28*)  ]/sin22\Jj 


3l/3a  =  0 


m3L/3B  =  cos2a /sin  (a+9*+i^) 

m3L/3a  =  -  [B  (^sin  2a  +  3y/3b)  +  H  3y/3h]  /  sin  (a+9*+ijj)  (B-4) 

m  3l/3h  =  (3Y/3H)  cos  (a+9*)  /cos^ 

m  3l/3’^  *  -B  (m  3L/3B)  /  tan  (a+0*+iii) 

+H(m3L/3H)  [  cosij;  +  ( 3Y/3H)  sin  (a+2t|>)  /2cosi^] 

The  dispersion  in  wavelength  is 

3Y/3A  =  -[b(3y/3b}2/  cos2  (Q*+\p) 

2  (B-5) 

+  H(3y/3h)  /  sin  2C'  sin{a+20‘)]  /2D  cos  9* 

B-l-b.  Special  case:  0  (Conturie  et  al  1981) . 

B-l-b-1.  Analytic  solution,  in  order  to  solve  Eq. (B-l)  one  needs  to 
have  at  least  as  many  known  lines  in  the  spectrum  as  there  are  unknown 
parameters.  In  the  general  case  it  is  best  done  on  the  computer. 

However  one  often  chooses  crystals  cut  along  the  direction  of  the 
diffracting  planes,  i.e.  ip  »0.  Eq. (B-l)  reduces  to  an  expression 
analytically  easy  to  handle: 

tan9  =  -tana  +  B/(Y-A)  (B-6) 

In  order  to  solve  for  the  three  parameters,  A,  B  and  tana  ,  one  needs  to 
know  three  lines  in  the  spectrum.  Calling  the  coordinates  of  these 
lines  (Y^,0^),  '  <^3 f 63) »  and  adopting  for  the  indices  the 

convention  that  i+3  stands  for  i,  we  have 


A  -  I  Y.[tane.<Y  -V  ))  /  I  t^.Y^-Y  ) 

1=1,3  1=1,3 


B  =  -  II  (Y±-Yi+1)  (tan9i-tan0i+1)  /  [  I  Y^+2  (tanS^-tanQ^^)  ] 
i=l,3  i=lr3 


tana  =  -  Z  tan9^ [Y^ (tan9^+^-tan0^+2) ]  /  Z  Y^  (tand^^-tanS^^) 
i=l,3  i*l,3 

The  partial  derivatives  are 

3A/3V.  -  (Yi+1-A>  l*i+2-A)  /  (Y^j-Y.)  <Vi+2-Vi) 

SB/SYj  /  ‘WV'VrV 


3tana/3Yi  -  -B/<Yi+1-Y.  >  (Y^-Y^ 


If  errors  AY^  (i=l,3)  are  made  in  the  true  positions  of  the  three 
reference  lines,  the  error  in  the  deduced  wavelength  of  an  unknown  line 
(Y,0 )  can  be  obtained  by  differentiating  Eq. (B— 6 )  with  respect  to  Y^ 
(i=l,3).  The  result  is 


A (tan0) 


-B  Z 
i-1,3 


(Y-Vi+2>  lVi+2-Ii*l)iYi 


/  (y-a)  n 
i=l,3 


B-l-b-2 .  Discussion  of  the  solution.  We  have  not  included  the 
2D-spacing  of  the  crystal  among  the  parameters.  Its  value  is  usually 
well  known  and,  in  addition,  the  wavelength  determination  is  insensitive 
to  an  error  in  it.  Differentiation  of  the  Bragg  tx.;:-.ition  and  of 
Eq. (B-6)  leads  to 

3A/3 (2D)  -  X/2D  +  2D  cos* 9  3tan0/3(2D) 


* 


3tan9/3  (2D)  -  -3  tana/  3  (2D)  +  {3b/3(2D)  +  [B/  (Y-A)  ]  3a/3  (2D)  }/ (Y-A) 
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For  the  three  reference  lines,  3A^/3(2D)  *  0  (i*l,3),  and  one  can  solve 
for  the  quantities  3A/3(2D),  3B/3(2D),  and  3tana/3(2D).  The  resulting 
error  in  the  measured  wavelength  is  given  by 

3A/3 (2D)  =  (2D)-3  n  (A-X. cose/cose. )  (B-12) 

i-1,3  1  1 

-4 

In  most  cases  the  right-hand  side  of  Eq. (B-12)  is  very  small  (£10  ), 
and  the  value  of  the  2D-spacing  need  not  be  known  with  great  accuracy. 

Since  the  parameter  tana  is  easy  to  measure  and  does  not  change 
with  the  position  of  the  spectrometer,  we  single  it  out  and  show  here 
that  the  measurement  of  tana  with  only  modest  precision  can  improve 
considerably  the  precision  in  determining  the  wavelength  of  an  unknown 
line.  Starting  with  Eq. (B-6)  and  assuming  that  tana  is  known,  we  solve 
for  A  and  B: 

A=  [Y2(tan62  +  tana)  -  Y1(tan91  +  tana)  ]  /  (tan©2  -  tan^) 

(B-13) 

B  =  ~  I  (Y2~Y^) / (tan02  -  tan9j, )  ]  (tan02  +  tana)  (tanS^  +  tana) 


Their  partial  derivatives  are  (with  the  convention  i+2*i) 

3a/3y.  *  (tan0.  +  tana)  /  (tan0.  -tan0.  ,) 
ii  i  l+l 

(B-14) 

3b/3y;.  =  B  /  (Yi-Y.+1) 


The  differentiation  with  respect  to  tana  of  Eq. (B-6) ,  where  A  and  B  are 
functions  of  tana  (see  Eq.(B-13)),  allows  us  to  compute  the  error  in 
the  determination  of  the  wavelength  of  an  unknown  line  (Y,0)  for  a  given 
error  in  tana  : 

2  3/2  2 

3X/3  tana  =  -2D(1-(X/2D)  ]  (Y-Y^)  (Y-Y2>/ (Y-A) 


(B-15) 


An  error  of  1°  in  a  will  cause  an  error  in  the  wavelength  AX/X  of  %2 

-4 

10  .  The  error  in  A  and  B  from  Eq. (B-14)  is  about  one  order  of 

magnitude  smaller  than  that  given  by  Eq. (B-8) .  The  error  in  the 

determined  wavelength  that  is  due  to  errors  in  (i*l,2)  is  obtained  by 

differentiating  Eq. (B-6)  and  by  using  Eq. (B-14) : 

2  3/2  2 

AX  =  2D[1-(X/2D)  ]  B[(Y-Y2)AY1-(Y-Y1)AY2]/[(Y1-Y2) (Y-A)  ]  (B-16) 


B-l-b-3.  Dispersion  characteristics.  When  Eq.  (B-5)  reduces  to 


3Y/3X  =  -Bco^a  /  2D  cos6  sin  (a+8) 


(B-17 


It  is  convenient  to  use  the  dimensionless  quantities  U«(Y-A)/B,  A*X/2D, 
E*l/A  to  compute  the  dispersions  in  wavelength  3U/3A  and  wavenumber 
3U/3E.  Results  are  summarized  in  Table  B-l  and  illustrated  in  Figures 
B-2-a  and  B-2-b.  The  flat-crystal  spectrograph  is  not  a  constant 
dispersion  device  (neither  in  wavelength  nor  in  frequency) ,  which 
precludes  extrapolation  far  outside  of  the  range  where  known  lines  lie. 
For  any  angle  0°<a<90°,  3U/3A  goes  through  an  extremum  for  some 
wavelength  A*Aq  (ex)  ,  while  this  happens  to  3U/3E  only  for  negative 
values  of  a.  3U/3A  varies  less  when  a  is  large,  3U/3E  varies  less  when 
a  is  small.  The  dispersion  in  wavelength  is  relatively  constant  around 
the  extremum,  where  a  linear  interpolation  might  yield  acceptable 
precision.  To  estimate  the  error  in  using  a  linear  interpolation  we 
develop  the  function  U(A)  in  a  series  expansion  around  the  extremum 
A  «Aq  up  to  third  order  (the  second  order  term,  of  course,  vanishes). 

We  finally  get  for  the  error 

(A-A*)/A  =  -(A-A  )  (A-A  )  (A+A  +A  -3A  )  .  (1-A2 /2)/l2A(l-A2  )2  )  (B-18; 

X  z  J.  2  O  O  O 


Dispersion  in  wavelength  (A  =  X/  2D) 

Dispersion  in  wavenumber  (Z  =  2D/X) 

3U/3A  =  -U-A2)_,|(A  +  /l-h2  tana)-2 

3U/3Z  =  E  (Z2-l)~’,(l  +  A*-l  tana)-2 

2  2 

3  U/3A  =  0  for  A=AQ  such  that 

32U/3E2  =  0  for  £=EQ  such  that 

tana  =  (1-3  Aq/2)  / (3  A  /2)  /l-A* 

tana  =  -[(2  + 1)  /l*-l  ] 

It  corresponds  to  a  minimum  of  -3U/3A  : 

It  corresponds  to  a  minimum  of  3U/3I  : 

[3U/3A]  .  =  -9A2/4/£aT 

[3u/3E’r  =  (2  Z2  +1)2/ 4  I3  /l2-l 

0 

'0 

Also  U(A  )  *  (3Aq/2)A^ 

Also  U(E  )  =  (1+  1/  2I2)  /l2-l 

3U/3A  •+•  when  A  -►  -sina 

3U/3E  -*■  +°°  when  I  ■+■  -1/sina 

3U/3A  -*■  when  A  -*■  1 

3U/3E  -*•  +°°  when  I  -*■  1 

3U/3A  -1/tan2  a  when  A  -*■  0 

3U/3E  -*■  0  when  E  +00 

(except  when  a=0.  Then  9U/3E  •*  1) 

Table  B-l  :  Dispersion  Characteristics  of 
a  Flat-Crystal  Spectrograph, 


Wavelength 

dispersion 
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Here  A^,  are  two  known  wavelengths,  A  is  the  actual  value  of  an 
* 

unknown  line  and  A  the  wavelength  of  the  same  line  obtained  by  linear 
interpolation  between  A  and  Eq. (B-18)  is  valid  in  a  limited  range 
around  A  ,  i.e.  in  the  region  where  the  curve  3U/3A  is  reasonably 
symmetric.  Outside  that  region  the  Taylor  expansion  must  be  pushed 
further.  The  upper  bound  of  j  A-A*[/A  for  any  wavelength  A£IA1,A2] 
within  the  range  of  validity  of  Eq. (B-18)  is 

I  A— A*  | /A  <  {  (1-A2o/2)/[4A(1-A2o)2]}  .Max  {  |A  -.  Aj3,  |A2-Aq  j3  }  (B-19) 

We  plot  in  Figure  B-2-c  for  each  angle  a  a  range  of  wavelengths  within 
which  linear  interpolation  will  yield  a  relative  error  smaller  than 


When  the  angle  a  is  known  accurately,  it  is  convenient  to  rewrite 
Eq. (B-6)  as  follows: 

'V/  a, 

Y  ■  A  +  Bcotania+9)  (B-20) 

<v,  %  2 

where  A*A  +  B.sina  .cosa  ,  B«B.cos  a  (see  Fig.B-l-b).  The  graph  Y  versus 
cotan (a+9)  is  a  straight  line  and  a  linear  interpolation  is  possible  for 
any  value  of  a  .  Also  the  dispersion  in  wavelength  9Y/3X  »-&/sin2 (a+6) 
derived  from  Eq. (B-20)  remains  valid  for  a  >90°  because  B  is  finite 
(while  B  becomes  infinite) . 


Curved-crvstal  s Dec t roqr a oh 


Mechanical  constraints  or  the  risk  of  damage  from  target  blow-off 
■ay  require  that  the  spectrograph  be  set  at  a  large  distance  from  the 
source,  reducing  the  range  of  angles,  hence  wavelengths,  incident  on  a 
flat  crystal.  The  spectral  range  may  be  increased  by  using  a  curved 
crystal.  The  relevant  geometry  is  presented  on  Figure  B-3.  We  limit 
our  analysis  to  the  case  where  the  diffracting  planes  are  parallel  to 
the  surface  of  the  crystal.  Five  independent  parameters  are  needed  to 
define  the  problem: 

CN  ■  r  :  radius  of  curvature  of  the  crystal  (center  C) , 

CS  ■  s  :  distance  from  C  to  the  source  S, 

CH  ■  h  j  distance  from  C  to  the  film  5^ 

(CH,CS)«  ot  :  (defines  the  orientation  of  the  film), 

OA  •  y  :  (defines  the  origin  0  on  the  y-axis  along  the  film) . 

Once  again  the  X-rays  are  diffracted  according  to  Bragg's  law, 

X  -2D.sin0  ,  where  D  is  the  interplanar  spacing  of  the  crystal.  One 
gets  from  Fig.B-4 


sinB  *  (r/s)  cos© 


y  *  y  +  [h  sin(8+20)/cosa  -  r  cos©]/ cos  (a+B+20) 


L  *  (2s  tan9  cos  (8+0)  +  (y-y)  cosa]  /sin(S+20) 


y  defines  the  position  of  the  line  of  wavelength  X  on  the  film.  L  is 
the  distance  (STS')  the  ray  of  wavelength  X  travels  from  the  source  to 
the  film.  If  a  densitometer  trace  of  magnification  m  is  used,  one  can 


make  the  substitutions:  Y«my,  R-mr ,  S*ms,  H»mh,  Y*my.  The  sensitivity 


of  Y  and  L  to  the  five  parameters  is  measured  by  the  partial 
derivatives: 


3Y/3R  =  -(S/R)  3Y/3S  +  (m  3L/3H)  cosQ 
3y/3s  =  -  (tanB  /S)  (3Y/3a  +  H  /  cos2  a) 

3Y/3H  =  -  (m  3L/3H)  sin(8+29)  /  cosa 
3Y/3a  =  (Y-Y)  tan  (a+8+20)  +  H  tana  3Y/3H 
3y/3y  =  1 

and 

m  3L/3R  =  {cosa  [Scos(a+20)  -  H]/cos2(a+8+26)  -  m  L  cos  (8+29) 
-  2S  tanQ  sin  2  ( 8+0 )  )  tanB  /  R  sin  ( 8+20 ) 

m  3l/3s  =  (L  -  R  3L/3R  -  H  3L/3H)  /  S 

m 3l/3h  =  -  1/  cos (a+8+20) 

m  3L/3a  ■  [R  cos9  -  H  sin(a+B+20)  3  /  cos2  (a<-3+20) 

3L/3Y  =  0 

The  dispersion  in  wavelength  is 


(B-24) 


(B-25) 


3Y/3X  *  (S  (tan0  -  2 /tanB)  3Y/3S  +  R  sin0/cos  (a+8+20)  ]  /  2D  cos0  (B-26) 


Figure  3-5  shows  the  geometrical  set-up  of  an  EUV  grating 
spectrograph.  The  grating  acts  as  both  the  dispersive  element  and  a 
1-dimensional  imaging  element.  For  each  wavelength  the  image  of  the 
entrance  slit  is  formed  on  the  film.  There  is  focussing  only  in  the 
direction  parallel  to  the  Rowland  circle's  plane.  One  gets  from  Fig. 3-5 

k  X/p  =  cosa  -  cost(y+a)/r]  =  cosa  -  cosB^  (B-27) 


L  *  r  (sina  +  sinB) 


(B-28) 


p  and  r  are  the  spatial  period  and  radius  of  curvature  of  the  grating. 
Positions  y  on  the  film  are  measured  with  respect  to  an  arbitrary  point 
0  whose  distance  to  the  center  of  the  grating  along  the  Rowland  circle 
is  a.  a  is  the  angle  of  incidence  and  6^  the  angle  of  diffraction  for 
the  wavelength  X  in  the  kth  order.  L  is  the  distance  (SCS')  the  ray  of 
wavelength  X  travels  from  the  entrance  slit  to  the  film.  If  a 
densitometer  trace  of  magnification  m  is  used,  one  can  make  the 
substitutions:  Y»my,  A*ma,  R»mr.  There  are  three  independent  parameters 
(A,  R,  a).  The  sensitivity  of  Y  and  L  to  these  parameters  is  measured 
by  the  partial  derivatives: 


3y/3a  -  -  1 

3y/3r  *  B  (B-29) 

3Y/3a  *  R  sina  /  sinB 


3L/3A  =  0 

m  3L/3R  =  sina  +  sinB  (B-30) 
m  3L/3a  =  R  cosa 


The  dispersion  in  wavelength  is  3Y/3X  =  kR/psinB 


B-3-b.  Discussion. 


In  order  to  solve  for  the  three  parameters  A,  R,  a,  one  needs  to 
know  three  lines  in  the  spectrum.  Let  us  call  the  coordinates  of  these 
lines  ^Y2'X2^  an<*  (*3*^3)  •  The  presence  of  trigonometric 

functions  in  Eq.B-27  does  not  make  possible  an  explicit  solution  for 
each  parameter/  similar  to  Eq.B-7.  Using  a  cosine  expansion,  one  gets 
the  following  results  in  first  order: 


=  -*5 


i=l»3 


ai+2“Xi+l)Yi  / 


i*l,3 


Ui+2-Xi+l)Yi 


=  ( -P/2 


i=L3 


(Yi+2-Yi+l)  / 


i=l»3 


(Y .  -Y .  . ) X . } ' 
1+2  l+l  1 


(B-31) 


I 

i=L3 


(Yi+2-Yi+l) 


'2  +  p  ot 


=  0 


However,  it  is  best  to  solve  for  an  exact  solution  on  the  computer. 

When  one  or  more  parameters  are  known  from  previous  recordings  with  the 

same  instrumental  set-up,  the  knowledge  of  less  than  three  lines  can  be 

enough.  The  algebra  to  compute  the  remaining  parameters  becomes 

straightforward.  It  is  not  presented  here. 

We  have  not  included  among  the  parameters  the  spatial  period  p  of 

the  grating  (p^)*ioVn;  n*number  of  grooves/mm).  The  expression  dVdn 

was  derived  for  each  possible  combination  of  free  parameters.  The 

results  are  summarized  in  Table  B-2  where  6^,  B ^  are  the  diffraction 

angle (s)  of  the  reference  line(s)  used  in  the  analytic  solution: 

8j»(Yj+A)/R.  The  table  gives  in  each  case  the  lower  order  term  (but  by 

far  the  dominant  one)  of  the  exact  solution.  Note  that  dVdn  goes  as 
2  2 

8  /2 1  (or  a  /2I)  except  in  cases  "2-3"  and  "1-2-3"  where  it  goes  as 


-[1-  (B1/3)2  ]  B2/2 


V  X3 


X3'  X1 


X.,  X_ 


Xl»  X2'  X3 


(i  +  z  e./B) [  n  (i - b./B)]  b  / 4 ! 

i=l,3  1  i=l,3  1 


X^  =  m  a  =  A 
X2  =  m  r  =  R 


X3  =  a 


S,  B^  in  radian 


Table  B-2  :  Effect  on  the  Line  Identification  of 

an  Error  in  the  Grating'  Spatial  Period 
in  a  EUV  Spectrograph. 
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3*/4i.  In  first  approximation  a  change  An  in  the  number  of  grooves  per 


mm  of  the  grating  does  not  affect  the  relationship  X«-+Y  if  it  can  be 
compensated  by  changes 

AA  =  0  ;  AR  =-  (R/2)  An/n  ;  Aa  =  [tan  (a/2)]  An/n 


This  of  course  is  possible  only  if  R  and  a  are  both  free  parameters. 

The  three  quantities  m,  r  and  n  are  often  reasonably  well  known  in 
advance  but  the  addition  of  just  one  free  parameter,  namely  R=mr,  takes 
care  of  the  uncertainty  on  all  three  at  the  same  time. 


B-4.  Least-squares  optimization. 


B-4-a.  Analysis. 

When  the  number  (N)  of  known  lines  in  the  spectrum  is  larger  than 

the  number  (M)  of  unknown  parameters,  one  can  improve  on  the  accuracy  of 

the  analytical  solution  given  above  by  doing  a  least-square  fit  with  the 

free  parameters  (j=l,M) .  The  errors  in  the  parameters  and  in  the 

wavelength  determination  are  now  determined  by  a  statistical  method. 

Let  us  label  the  N  reference  lines  as  follows:  Y.±AY.,  X.+AX.  (i=l,N) , 

1  l  i  l 

where  A^  corresponds  to  the  uncertainty  in  the  position  of  the  line 
center  for  the  ith  line  and  AX^  refers  to  the  accuracy  of  the  tabulated 
wavelengths.  We  define  the  real  position  Y^+6Y^  corresponding  to  the 
wavelength  X^,  by  a  normalized  probability  distribution 


n  *  6Yi 

P /I  <$v  \  *  -  exp  -  - - - 

i'  i  AYA  +  (3Y/9X)iAXi  [_AYi  +  (3Y/9X)iAXi 


(B-32) 


We  further  assume  that  the  errors  AY^'s  made  in  the  measured  position  of 


\-:y% 


the  N  lines  are  independent  of  each  other,  i.e.  <  6^.  6Yif>  °* 

Using  Eq.B-32,  we  have 


(B-33) 


The  knowledge  of  j=l,M,  and  j'*l,M)  allows  us  to  find 

-  the  position  of  a  given  wavelength  f(Xj),X], 

-  the  wavelength  of  an  unknown  line  X=f  1 E(X^ ) , Y ] ,  (B-34) 

-  in  each  case,  the  distance  L»g[(Xj),X], 

source-film  (crystal  spectrograph)  or 

entrance  slit-film  (grating  spectrograph) 
and  to  compute  the  rms  errors  associated  with  the  solution  of  Eq.B-34: 


rms (5x J 

■  ''V1 

a 

rms  (6y) 

n 

ii 

^  M 

se 

Z  (3y/3x.)  (3y/9x.,)  S. } 
j'=l,M  3  3  33 

b 

rms (6X) 

-  (9X/9Y) 

rms (5Y) 

c 

rms  (6D 

=  z 

j=l,M 

Z  (3l/9x.)  (9l/9x.,)  S.  .,  } 

j'=lM  3  3  33 

d 

The  partial  derivatives  9Y/9X^,  9L/9X^  (j*l,M)  needed  to  compute 
quantities  B-35-b  and  B-35-d  were  given  in  the  previous  sections:  see 
Eq.B-3  and  B-4,  Eq.B-24  and  B-25,  Eq.B-29  and  B-30. 

FORTRAN  programs  were  written  to  carry  out  these  calculations. 

Their  filenames  are  CRYSPEC  (for  the  flat-crystal  spectrograph) ,  CURVEX 
(curved-crystal)  and  GRASPEC  (grating).  The  parameters  derived  from  the 
analytical  solution  with  M  known  lines  are  used  as  initial  values  for 
the  least-squares  fitting.  Let  us  define  F^  *  Y^  -  f [ (X ^ ) , X^ ]  . 


The  program  searches  for  the  values  of  X^  (j=l,M)  that  minimize  the 

2 

quantity  E  (i*l,N).  It  computes  the  matrix  ||  9X^/3 Y^||  whose 
components  are  needed  in  Eq.B-33: 

H3X./3Y.H  -  (J^)^  jJ.  (B-36) 

where  J  is  the  jacobian  matrix  ||  9F^/3 x ^ 1 1  and  Jfc  is  the  transposed 
matrix.  The  program  output  includes  the  solution  to  Eq.B-34  and  B-35 
and  the  value  of  the  local  dispersion  (3Y/3A)  for  each  line. 

B-4-b.  Examples. 

The  program  CRYSPEC  has  been  tested  with  a  large  number  of  known 
spectra  in  the  region  1-20  ft.  The  uncertainties  AY^  and  AA ^  in  the 
reference  lines  were  set  equal  to  the  half-width  at  half-maximum  and  to 
one  unit  of  the  last  digit  listed  in  Kelly  &  Palumbo  (1973), 
respectively.  The  difference  between  the  wavelength  determination  from 
the  spectrum  with  our  method  and  the  tabulated  value  was  always  within 
the  error  bar  (AA  /\<,  10  predicted  by  the  program.  An  example  is 
shown  in  Figure  B-4.  A  Chlorine  spectrum  was  recorded  with  a  P.E.T. 
crystal  (tetragonal:  a^»a2*6.16  ft;  a^*8.74  ft).  The  spectra  produced  by 
two  sets  of  diffracting  planes  were  recorded  simultaneously.  One  set 
was  (0,0,2),  parallel  to  the  surface,  2D»8.74  ft.  It  allowed  the 
determination  of  the  parameters  A,  B  and  a  (see  Sec.B-1).  The  remaining 
lines  were  then  found  to  have  been  diffracted  by  the  set  (0,1,3), 
ip  m2§  20',  2D*5.27  ft.  The  parameter  H  was  also  computed. 


The  program  GRASPEC  has  been  used  to  identify  many  spectra  in  the 
range  45  ft-200  ft,  from  Beryllium  (Z»4)  to  Iron  (Z*26) ,  with  gratings  of 


Design  and  fabrication  of  the  cylindrical  optics 


The  beam  produced  by  the  Glass  Development  Laser  and  the  lens  used 
to  focus  it  on  target  are  well  characterized  in  both  the  IR  and  the  UV 
configurations  (Seka  et  al  1980  &  1981).  In  some  of  our  experiments  we 
needed  to  add  a  blast-shield  to  protect  the  lens  from  target  blow-off 
and/or  a  cylindrical  optics  to  form  a  line-focus.  This  appendix 
describes  a  ray-tracing  code  written  to  evaluate  how  the  blast-shield 
affects  the  intensity  distribution  on  target,  in  a  point-focus  or  a 
line-focus  geometry.  It  also  describes  several  options  for  the 
cylindrical  optics,  the  design  that  was  chosen  (a  pair  of  weakly 
diverging  corrector  plates)  and  the  manufacturing  process. 

C-l.  Ray-tracing  code. 

In  high  intensity  experiments  on  flat  targets,  the  plasma  blow-off 
occurs  mainly  in  the  direction  of  the  incident  laser.  When  a  cooling 
plate  is  added  in  front  of  the  target,  the  expanding  plasma  can  produce 
hot  debris  which  are  then  thrown  toward  the  focusing  optics,  damaging 
its  coating.  An  inexpensive  glass  plate,  inserted  at  an  angle  to  avoid 
back  reflection  and  ghost  formation  along  the  axis  of  the  laser,  is  used 
as  a  shield  (see  Figure  C-l-a) .  The  code  (filename  LINFOC)  can 
propagate  a  uniform  parallel  beam  through  a  2-focus  lens  followed  by  a 
glass  plate  (i)  of  any  index  and  thickness,  (ii)  with  any  tilt  in  any 
direction.  Since  the  introduction  of  the  glass  plate  and/or  the 
astigmatism  of  the  cylindrical  optics  are  the  main  factors  influencing 
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Figure  C-l  :  Optics  to  Produce  a  Line-Focus. 

A-  TILTED  BLAST-SHIELD, 

b-  (geometry  for  the  code  LINFOC) 

c-  CYLINDRICAL  CORRECTOR  PLATES. 


s  *  cosu 
z 


The  tilt  of  the  glass  plate,  assumed  plane  and  parallel,  is  defined  by 
the  angle  Y  between  its  normal  and  the  optical  axis  Z,  while  its 


For  an  ASE  device  in  the  soft  X-ray  region,  one  needs  to  produce  a 

line-shaped  plasma  of  aspect  ratio  p£lO:l.  A  smaller  ratio  would  make 

it  difficult  to  observe  any  angular  dependence  because  plasma  expansion 

tends  to  blur  the  geometry  of  the  original  focus.  It  is  not  easy,  at 

least  in  our  experimental  set-up,  to  create  a  plasma  of  transverse 

dimension  smaller  than  80  ym.  The  minimum  length  of  the  line  focus 

should  therefore  be  ^  1  non.  The  maximum  length  depends  on  the  available 

laser  energy  and  the  atomic  number  of  the  target  material.  A  population 

inversion  in  the  Helium-like  stage  of  Aluminum  requires  about  2xl014 

W/cm^  of  UV  light  or  4xl014  W/cra^  of  IR  light  on  target.  A  1  mm  x  100 

ym  focal  spot  would  require  100  J  of  UV  light  or  200  J  of  IR  light  in 

500  psec.  It  is  close  to  the  limit  capability  of  the  GDL  laser  when 

upgraded  with  active  mirrors  (see  Chapter  3) .  Targets  of  lower  atomic 

13  2 

number  such  as  Fluorine  require  intensities  in  the  10  W/cm  range, 
making  it  possible  to  extend  the  length  of  the  line  to  several  nm. 

There  are  a  number  of  designs  to  form  such  a  line-focus.  Limiting 
ourselves  to  a  maximum  of  three  optical  elements,  to  reduce  reflection 
losses  and  ghost  problems,  we  can  consider: 

a)  one  single  lens  with  two  different  foci, 

b)  two  strong  cylinders,  in  translation  or  rotation, 

c)  the  main  aspheric  lens  (the  lens  used  in  point  focus  experiments) 
followed  by  a  thick  blast-shield  with  a  large  tilt  angle. 


d)  two  weak  rotating  cylinders,  before  the  aspheric  lens. 


The  first  two  options  are  difficult  to  manufacture  in  large  size. 
Also  in  option  a  the  length  of  the  line  cannot  be  tuned.  We  did, 
however,  use  an  existing  lens  of  that  type  in  preliminary  experiments 
with  IR  light  (see  Chapter  4) . 

Option  c  takes  advantage  of  the  strong  astigmatism  introduced  by 
the  blast  shield.  Figure  C-2,  produced  by  the  code  LINFOC,  illustrates 
how  one  can  form  a  line-focus  of  aspect  ratio  =10:1  using  that  method. 
The  length  of  the  line  can  be  tuned,  though  in  a  limited  range,  either 
by  changing  the  tilt  angle  y  (which  does  not  require  another  plate)  or 
the  thickness  t.  There  is  a  large  lateral  shift  (6  ■  several 
centimeters)  of  the  laser  axis.  One  may  expect  the  non-uniformity  of 
the  geometrical  image  to  be  smoothed  by  diffraction  in  the  beam  or  the 
optics,  though  it  is  not  a  requirement:  the  focus  is  dissymmetric  enough 
for  our  purpose  and  its  central  gap  may  actually  help  the  cooling  of  the 
plasma  (Silfvast,  et  al,  1979-a) . 

Option  d  is  the  most  attractive.  It  is  easier  to  manufacture  than 
options  a  and  b  and  more  tunable,  less  aberrated,  and  less  massive  than 
option  c,  with  no  beam  walk-off.  It  is  described  in  the  next  section. 

C-3  Design  and  manufacture  of  a  cylindrical  lens. 

Figure  C-l-c  shows  the  combination  of  lenses  chosen  to  produce  a 
variable  line  focus.  S  is  the  aspheric  lens  used  in  point-focus 
experiments  (focal  length  f),  and  are  two  identical  cylindrical 

lenses  with  very  weak  power  (l/|f  |  ■  1/ 1 f  |  *  1/f  «  1/f ) •  Neglecting 

Cj  c2  c 

the  spacings  between  the  lenses,  the  length  of  the  line  focus  is 
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LASER  INTENSITY  DISTRIBUTION 
ON  A  FLAT  TARGET 
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Figure  C-2. 


L(9)  =  2*  (f/f  )  cos  29  (c-2) 

c 

where  D  is  the  diameter  of  the  incident  beam  and  9  the  half-angle 

between  the  axis  of  the  two  cylinders.  The  relationship  between  f  and 

c 

the  cylindrical  sag  at  the  center  of  one  of  the  lenses  is  as  follows: 

-l/fc  *  l/fc  =  (n-l)/R  (C-3) 

sag  =  R(l-cosa)  -  R  ot2  /2  =  $/8(n-l)f  (C-4) 

c  c 

where  R  is  the  radius  of  curvature  of  the  cylinder.  Therefore 

sagc  *  Lmax/16(n“1)  f#  (c_5) 

Since  the  absolute  value  of  the  overall  focal  length  (F=  f)  of  the 
three-lens  optics  is  unimportant,  the  cylindrical  lenses  need  not  be 

perfect.  If  C.  and  C,  have  some  spherical  power  (1/f  ,  l/f_  «  1/f) , 

A  S2 

the  focal  length  is  changed  by  AF  a (1/f  +  1/f  )f2.  The  power  1/f 

Sj  S2 

results  in  a  spherical  sag  in  the  ith  lens: 

(A  sag)  =  <t>/8(n-l)f  *  (AF)  $ /8  (n-1)  f*  *  (AF)  .  /8  (n-1)  f  #  (c-6) 

i  i  1  J- 

If  one  tolerates  an  overall  uncertainty  ±(AF)  on  F,  the  tolerance  on 

max 

the  spherical  sag  for  each  "cylindrical  lens”  is 
(Asag)  -  (AF)  /[(n-l)(4f#n  (07) 

S  max 

To  have  the  flexibility  of  working  with  low-z  targets,  we  set  the 
tunability  requirement  to  be  L  ■  0  to  4  mm  (L  -  2  mm  when  only  one 
cylinder  is  used).  With  n*1.5  and  (AF)max=l  cm,  the  specifications  on 


each  lens  ace  as  follows: 


f# 

sagc 

(  A  sag) 

s 

3 

166  ym 

139  ym 

6 

83  ym 

35  ym 

12 

42  ym 

9  ym 

The  UV  beam  and  focusing  optics  available  to  us  were  such  that  D  *  13  cm 

and  f -number*  12,  hence  the  idea  arose  of  achieving  the  small 

cylindrical  sag  required  by  bending  and  polishing  a  flat  plate  of  fused 

silica.  The  process  is  described  on  Figure  C-3.  The  dimensions  of  the 

plate,  not  including  the  edges,  were  5.7"x5.7"x0.3" .  The  bending  knife, 

which  ran  the  length  of  the  plate,  was  raised  by  the  motion  of  two 

sliding  wedges.  Each  step  was  controlled  with  a  Fizeau  interferometer. 

The  main  difficulties  of  that  process  are 

(i)  to  polish  flat  a  large  square-shaped  surface  (step  1:  lower 

surface,  step  3:  upper  surface) .  It  is  difficult  to  maintain  the  same 

rate  of  wear  at  the  corners  of  the  plate  and  at  the  center.  Fortunately 

the  knife  runs  along  a  transverse,  not  a  diagonal  direction  and  in  the 

end  the  lens  is  used  in  a  circularly  symmetric  beam.  The  quality  of  the 

corners  affects  neither  the  symmetry  of  the  bending  nor  the  quality  of 

the  line-focus.  Along  the  direction  of  the  knife  the  flatness  of  the 

plate  must  be  within  about  a  wave,  the  mechanical  quality  of  the  knife 

becoming  the  limiting  factor.  The  overall  polish  must  be  smooth  enough 

for  the  deposition  of  anti-reflection  coatings  capable  of  handling  a 

high-power  beam.  In  our  experiments  the  power  load  at  0.351  urn  can  be 
2  2 

as  high  as  1  GW/cm  ,  or  0.5  J/cm  in  500  psec. 


MANUFACTURE  OF  A  CYLINDRICAL  CORRECTOR 
FOR  LINE-FOCUS  EXPERIMENTS 


1. 


The  lower  surface  has  been 
polished  flat.  The  plate 
is  Inserted  In  the  mount. 


2. 


The  edges  are  clamped. 
The  central  knife  Is  raised 
to  produce  the  appropriate 

aag. 


After  proper  blocking,  the 
upper  surface  is  polished  flat. 


The  plate  Is  unblocked. 
The  knife  Is  lowered, 
producing  a  concave 
(negative)  cylindrical  lens. 


Figure  C-3. 


(ii)  to  measure  the  sag  introduced  in  step  2,  which  determines  the 


length  of  the  line-focus.  The  bending  produces  a  saddle-shaped 
interferometric  fringe  pattern  because  initially  the  plate  is  slightly 
concave  at  the  center  (if  that  residue  was  convex,  the  bending  would  be 
less  even  and  the  pattern  more  distorted).  The  fringes  are  very  dense 
at  the  edges  (see  Eq.C-4).  Only  at  the  center  of  the  plate  can  their 
spacing  be  resolved  accurately  and  used  to  measure  the  sag  from  an 
extrapolation  based  on  the  assumption  of  a  circular  bend.  Good 
agreement  was  found  between  the  cylindrical  sag  set  inter ferometrically 
in  the  optical  shop  and  the  length  of  the  line-focus  measured  with  the 
alignment  beam  of  the  GDL  laser. 
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